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Some Properties of a Glass Used in Paper Manufacture 
Martin J. O'Leary and Donald Hubbard 


The chemical durability, heterogeneous equilibria, hygroscopicity, pH response, dis- 
persal, and thermal expansion characteristics of a fiber glass used in paper manufacture were 
investigated. The chemical durability of this glass was found to be very poor between pH 
2 and pH—2, apparently due to the high alumina and low silica composition. Above and 
below this range the chemical durability greatly improved, showing with the interferometer 
no detectable attack for exposures of 6 hours at 80° C either in high concentrations of H,SO,, 
or between pH 4 and pH 8. At alkalinities above pH 8 the durability again decreased as is 
characteristic of many silicate glasses. In addition, the heterogeneous equilibria at the 
solution-glass interface, using Ag(NH;).* ions as indicator, showed the fiber glass when 
leached at pH 4.1 to have an uneven distribution of migratable ions between the negatively 
charged glass surface and ambient aqueous solutions, to a greater extent than most commer- 
cial glasses. In accord with the very low hygroscopicity and poor chemical durability, 
electrodes prepared from this glass gave no definite response to hydrogen ion activity of the 
ambient solutions. The dispersal characteristics of the fibers in aqueous solutions follow 
the “critical” {points indicated by the chemical durability curve covering the range from 
pH—S8.1ltopH 11.8. Theannealing and softening temperatures as obtained from the respec- 
tive expansion curves were found to be higher for the fiber glass than for either Pyrex 7740 or 
Kimble N 51-A. 

The chemical durability data, the acid titration of sodium silicate, and the dispersal-pH 
series are consistent with the findings of O’Leary et al. in emphasizing pH 2.9 as a critical 
condition for the most successful production of paper from fibers of this glass. 


1. Introduction 


During the developmental stages in the production 
of paper from glass fibers [1,2]! it became evident 
that optimum conditions of strength and quality of 


durability, heterogeneous equilibria, hydroscopicity, 
pH response, dispersal, and thermal expansion 
characteristics of a fiber glass (E type, table 1) used 
in the production of paper; a general comparison 
was made with a few of the more familiar commercial, 


the finished product were obtained upon adjusting 
the hydrogen ion concentration of the suspended 
fibers to near pu 2.9 prior to the fabrication of the 


experimental, optical, and natural glasses. 





paper (ah the ydeoge an activity WHS | anus 1. Composition of ames eranged in he de fd 

strength of the finished paper was greatly ceduced. creasing percentage of SiO; 

Of considerable significance from an industrial 

standpoint was the observation that the fibers dis- o te 

persed rapidly at this optimum condition and dis- " «|2] |e Be &. | | 

tributed themselves evenly throughout the liquid zi. lel 2% (elisls|s] abla! 

medium, insuring uniformity of the resulting sheet. z Tie is) Slsizlalisis lola 

At high pH values the fibers showed a strong tend- 5 S/S ISI 2 (S188 l2 lais 

eney to remain in clumps and required a much longer - : A tad Desa 

time to disperse adequately for paper making [1,3]. so : oan See. | 0.1 
Some of the improved strength of the paper was B01. || 29 (tae) Rel ee | eee 88 

attributed to the cementing together of the dried AsO 0.005 0.03 0.5 ...--| 0.2] 0.5 

fibers by the swollen gelatinous silica-rich layer [3] piv : 

that is known to be formed on many types of silicate ~ “"h salea 14. 5 _ ms 

glasses during leaching in acid buffers [4 to 9]. Fur- ee n=-t- i 0.6 _— 

ther partial rationalization of the cause for this Cad | 0.9| ‘860 16.0 33.5 

increased strength at or near pH 2.9 was postulated Bao “s ei 

from the titration curve of sodium silicate by HC! wd 0.1 i. 7.2). 

[10], typical voltage departure curves of the glass HO ' 6 - waste 

electrode [11, 12], and the repression of swelling of — S/S G4) 506 6 oD 0.4) 3.6 

electrode glasses in the “super acid”? region below Ks0. 02/02 05) 40.._j1L8 7.7| 5.2 


pu 2 (5, 13). as 

In order to obtain a better understanding of the 
fabrication of paper from glass fibers, the present 
investigation was undertaken to study the chemical 


* Analyses by F. W. Glaze, J. Research NBS 26, 538 (1941) RP1394. 
l'ypical composition of obsidian, George W. Morey {30}. 
M. Dole {11}. i 

! Calculated from bath composition. 

e Approximate composition furnished by the manufacturer. 


‘ Typical composition, R. H. Bogue, The chemistry of portland cement (Rein- 
Figures in brackets indicate the literature references at the end ot this paper. | hold Publishing Corp., New York, N. Y., 1947). 








2. Experimental Procedure 


The determinations of chemical durability were 
made by the interferometer procedure which indi- 
cates in terms of fractions of a wavelength of light 
the surface alteration of the specimen after exposure 
to solutions under controlled conditions of time, 
temperature, and pH [9]. The light source used as a 
standard of length was the 5876 A line of the helium 
spectrum. For the range pH 1.8 to 11.8 the Britton- 
Robinson universal buffers were used [10, table 57c]; 
for the more acid exposures aqueous solutions of 
H.SO, were employed. All exposures were made at 
80° C and a durability comparison was made on a 
6 hr basis, although some of the exposures were 
limited to shorter periods. The surface alterations 
were plotted as positive values in the case of attack 
and negative values in the case of swelling. The 
specimens used in these experiments were polished 
cross sections of the various glasses whose surfaces 
were ground sufficiently flat to show reasonably 
straight interference bands when illuminated under 
a fused silica optical flat. In the case of the fiber 
glass the interferometer test specimens were prepared 
from ‘‘marbles’” of the glass used in fiber manu- 
facture. 

In the experiments on heterogeneous equilibria, 
the uneven distribution of the migratable ions, 
Ag(NH;,).* and Br-, at the solution-glass interface 
was demonstrated by potentiometric titration of the 
Agt and Br~ ions. The Ag:AgBr and glass electrodes 
were used as indicator and reference electrode, 
respectively [13, 14]. 

The hygroscopicity values represent the amounts 
of water sorbed by powdered samples, (approxi- 
mately 1.5 g that passed through a Tyler standard 
150 mesh sieve, equivalent to a U.S. Standard Sieve 
No. 140) upon exposure for 1 and 2 hr to the high 
relative humidity, approximately 98 percent, main- 
tained by a saturated solution of CaSO,2H,O at 
room temperature (6). 

The expansion measurements were made, using 
Saunders’ modification of the Fizeau interferometric 
method [15]. 

Glasses of the composition given in table 1 are 
those compared in this investigation. 


3. Results and Discussions 


3.1. Chemical Durability 


In figure 1, plotted from the data given in table 2, 


the chemical durability of the fiber glass is compared 
with the chemical durabilities of a Na,O-PbO-SiO, 
glass and with three typical optical glasses--BaC 
BSC 517, and F 620— over a range pH 2 to pH 
11.8 for exposure periods of 6 hr at 80° C. These 
four types are examples of glasses having poor 
chemical durability. Figure 2 and table 3 offer a 


O42, 


similar comparison with Corning O15 and_ three 
glasses of known chemical durability, Pyrex 7740, 
fused silica, and a natural glass (obsidian from 


Yellowstone Park). In these figures the fiber glass 
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FiGguRE 1. Chemical durability of the fiber glass used in the 


manufacture of paper compared withfthree typical optical 
glasse s, BSC 517, BaC 472, and F 620: and with an er peri- 
mental Na,O-PbO-SiO, (soda flint) glass. 


Data by Bernard R, Nebel, M. D, 


exhibits poor durability in the acid range below pH 4. 
However, it shows none of the conspicuous break- 
down of chemical durability between pH 5 and pH 9, 
characteristic of the optical glasses, particularly 
BaC 572. In the alkaline range the chemical dur- 
ability of the fiber glass is inferior beyond approxi- 
mately pH 8, as is characteristic of many silicate 
glasses. 

The decrease in the chemical durability indicated 
for the fiber glass below pH 4.1 made it desirable to 
obtain information at still lower pH values. The 
data in table 4 and figure 3 extend the curve for 
chemical durability in the acid range to pH —8.1 
by using aqueous solutions of H.SO,. Although it 
is obvious that such values as pH —8.1 cannot be 
interpreted in the Arrhenius sense [16], i. e., cannot 
represent concentration of hydrogen ions in moles 
per liter, nevertheless, the pH values calculated from 
measurements by the hydrogen-saturated KCl- 
calomel cell are plotted in order to extend the curve 
for convenience in representing the data. 
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TABLE 2. 


and an experimental Na,Q—PbO-SiQ, glass 
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Chemical durability of fiber glass used in the manufacture of glass paper compared with three typacal optical glasses 
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0 ———— oO of the resulting curve it is desirable to observe the 
appearance of the specimens (surface colorations) / 


after exposure at each pH (table 4 and fig. 3). From s 
9 pH 11.8 to pH 4, and from pH —2.4 to pH —8.1 thes on t 
surfaces are clear. The specimens exposed within — oper 
the range from pH 3 to pH 1.5 show interference the 
B colors (reflection-reducing films) on the surface, even tion- 


for those specimens at pH 0.75 and pH O that were migr 
etched and pitted. Such colors are normally the! the 














7 result of preferential leaching of cations from the — char 
glass, thus leaving behind a swollen silica-rich layer, of ¢ 
| similar to that exhibited by Corning O15 glass. and migr 
6 + by the optical classes under certain conditions and ¢ 
w . 
> (table 2) solut 
c one . > . 7 
e | Phe pu at which color first becomes evident for exces 
= St GR" GS CONEN Cnnuee the fiber glass coincides with the optimum pH 2.9; _ rise t 
° i . . . . . oon 
5s | for the production of paper from this glass. The — the s 
® fact that the interferometer disclosed attack instead mech 
| > . . , ~ . 
<a 4} of swelling for the fiber glass is the result of the of m 
S| partial sloughing off of the relatively weak silica-rich —_relati 
c r . P rie a‘ . . 
> | . = 5 laver produced initially [133]. The high ALO, high as We 
| °o ®@ ’ ° . . °° ML . 5 . 
” 3f ® a =< alkaline earth, low SiO, composition of the fiber glass inter’ 
| .o) . . . . . 
2 © | (table 1) is the main contributing factor in this and | 
m6 2:2 Corning breakdown This breakdown does not take place uselu 
2} =\z ala j o's in the presence of equal amounts of ALO; and high meth 
percentages of SiQ,, as exemplified by the behavior! chose 
of obsidian (fig. 2). That the effect is not centered ready 
J / around the CaO has been demonstrated in published pH, : 
8 | JA, + work on Na,O-CaQO-SiO, glasses [6]. This point is Fig 
a / further emphasized in table 5 and figure 4 in which the fi 
Oo as © + + found = —_ . y . 2 
ee at — is shown the weakening of chemical durability con- nugra 
s tributed by the presence of relatively large percent- the ¢ 
: | ages of ALO, in a glassy, high-alumina cement! (740, 
on ne ae ae ion ”. lk 
“te r 4 —<; 8 0 3 ia Clinker. ‘The data indicating the pH values at which glass 
the transition takes place for the solution and pre. — hydr« 
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The specimens showed very vigorous attack | swelling at pH 10.2 shown as a broken line. Very 2 + 
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3.2. Heterogeneous Equilibria 


Some insight into the nature of the ionic charge 
on the vlass surface and the disintegrating forces 
operating upon glass can be obtained by studying 
the heterogeneous equilibria appearing at the solu- 
tion-glass interface. The uneven distribution of 
migratable ions at this interface brought about by 
the presence of the nonmigratable, negatively 
charged silicate ions left after preferential leaching 
of cations must necessarily result: in an 
migratable positive ions appearing in the glass surface 
and of an equal excess of negative ions in the ambient 
solution The over-all result at equilibrium is an 
excess of electrolytes in the glass surface, which give 
rise to Osmotic pressure across the interface, causing 
the surface to swell as long as it Is chemically and 
mechanically durable [13]. The uneven distribution 
of migratable ions gives a ready indication of the 
relative chemical durabilities at pH 4.1 of the glasses, 
as well as a partial picture of the mechanism of the 
interfacial voltage potential between the solution 
and the glass surface. This procedure is especially 
useful under conditions in which the interferometer 
method is insensitive. The value pH 4.1 was 
chosen because the interferometer does not make a 
ready differentiation of many silicate glasses at this 
pH, as can be seen from figures 1 and 2. 

Figure 5 plotted from data in table 6 shows tha 
the fiber glass has a greater uneven distribution of 
migratable ions than the natural glass, obsidian, or 
the chemical and industrial glasses, such as Pyrex 
7740, fused silica, or Vycor. Corning 015, a special 
class whose voltages serve as an indicator of the 
hvdrogen-ion activity of aqueous solutions [11], also 
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TABLE 6. Uneven distribution of migratable Ag(NHe2)} tions, 
between the glass surface and the ambient solution for a fiber 
glass compared with obsidian and four special commercial 


glasses 
Excess g(NH ; ions in surface of 
specimens leached at pH 4.1, 80° C for— 
Glass 
0 hr 6 hr 24 hr 
M-eq/¢ M-eq/g M-eq/9 
V veor 0. 110-4 0.3x10-4 0. 7510-4 
Fused SiO l 1.4 3, 
Obsidian l 0.6 2.8 
Pyrex 7740 2 4.2 5.4 
Fiber glass 2 7.1 18. 2 
Corning 015 5.8 23. 0 30. 6 


has a larger uneven distribution of migratable ions 
than the other glasses. 

It is of interest to note the curves for obsidian and 
fused silica cross, as observed in earlier experiments 
114], showing that the rate at which negatively 
charged silicate ions are made available during leach- 
ing changes for these two glasses. Another interest- 
ing observation is that under the conditions of the 
test the sample of Vycor was superior to the specimen 
of fused SiQy,. 


3.3. Hygroscopicity 


Values of hygroscopicity (water sorbed from the 
vapor state) for the fiber glass compared with other 
glasses are given in table 7. Low values of hy- 
groscopicity are indicative of the ability of a glass 
to maintain a clear polished surface upon exposure 
to water vapor, a property different from chemical 
durability. The latter term is generally applied to 
the capacity of a glass to withstand attack by solu- 
tions [18]. The chemical durability of the fiber glass 
is inferior, although its hygroscopicity compares 
favorably with that of fused silica and Vycor. This 
low hygroscopicity is in accord with the high elec- 
trical resistance of the surface shown by the fiber 
glass at high relative humidities [19, p. 407]. 

TARLE 7 Chemical durability, heterogeneous equilibria, hygro- 
scopicity and pH response of fiber glass for comparison with 

a series of NagO-CaO-SiO» glasses and a few commercial and 


natural glasse 8 


Water sorbed 


Glasses Leached 6hr, 80° C at pH 4.1 
DH re- 
Surface Excess > sponse 
NiwwO | CaO | SIO alteration Ag(NH thr 2hr 
Fringes M-eq/a mg/cm? mg/cm? me/pH 
Ww a1 70 NDs $.1xl0-* 14 21 0 
15 70 NbD 5.3 20 $6 $1 
20) iT 70 ND 7.2 35 SS 59 
2 F 70 1/2+58 42 70 160 59 
2t t 70 28 i 80 185 59 
27 ; 70 6s 109 v2 219 57 
Ds 2 70 30S On) 108 249 55 
m9 ( *() 75S snot 135 314 42 
0 0 70 1268 <2t 152 379 19% 
Fused SiO: ND 1.4 6 10.5 0 
V veor NbD 0.3 8 11 0 
Fiber glass ND 7.1 9 11.2 0 
Obsidian ND 0. ¢ 8.9 3.9 0 
Australite ND l.¢ 10.1 15. 3 0 
Pyrex ND 4.2 18.5 26 : 
Corning 2/10-—S 23 61 125 9 


Symbols explained in table 2 
Values indicate that silica rich layers had sloughed off [13] 
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3.4. pH Response 


Glass electrodes prepared from samples of the 
fiber glass showed ho definite voltage response to 
hvdrogen-ion activity of buffer 
(table 7 This is in accord with previous observa- 
tions that the voltage of electrodes prepared from 
glasses of low hygroscopicity have always failed to 
satisfactory indicators of 
activity of aqueous solutions |[6. 11] 

An over-all picture of the position occupied by the 
fiber glass with respect to other glasses is given in 
table 7 and figure 6. In this figure the properties of 
chemical durability, uneven distribution of Ag 


(NH ‘ 


aqueous solutions 


serve as hvdrogen-ion 


» lons, and pH response are plotted as vari- 
ables dependent upon the hvgroscopicity as the 
abscissa. A Na0-CaO-SiO con- 
taining 70 percent of SiO, is included for comparison 
For qualitative orientation this series of reference 
glasses is divided into three utility groups based on 


SeTICS of 


classes 


commercial experience: superior, intermediate, and, 
inferior. It can be seen that the value obtained for 
the uneven distribution of Ag (NH ions for the 
fiber glass does not coincide with the trend Indicated / 
by the other classes This illustrates the fact that 
the chemical durability of this glass as determined 
from heterogeneous equilibria data is not as eood as 
the hvgroseopicity might imply 


3.5. Dispersal Characteristics 


The tendeney of the glass fibers to disperse and dis- 
tribute themselves uniformly throughout the aqueous 
medium ts illustrated in figure 7. Each member of 
the series consists of 0.25 ¢ of fiber immersed in 50 
ml of buffer at room temperature 
prepared for direct observation, and comparison with 
the chemical durability curve for fiber glass shown in, 

>. The solutions range from pH 8.1 to 


The series was 


figure 3. 
pH 11.8 obtained with H.SO, and Britton-Robinson 
universal buffer mixtures. 

It is important to observe the correspondence 
between the pu values indicated by the peptization 
series and some of the critical points shown by the 
chemical durability curve. Perhaps the most inter- 
esting observation is that wherever the glass surface 
after exposure remained clear, the fibers persisted in 
clumping. This is observed for solutions more 
alkaline than pH 4.1 and those more acid than pH 

1.49. In the acid range where interference colors 
appear on the surface of the exposed specimen, there 
is rapid dispersal and swelling of the fibers with the 
formation of gelatinous masses in the beakers for the 
region where the attack is the severest The practi- 
cal optimum found for the production of paper wit! 
this glass comes in the narrow range near pH 2.9 at 
which the first appearance of color begins and befor 
the region of severe acid attack starts 

Another interesting feature is that in fon? from 
pH 2.9 to pH 4.1 or higher and return, the dispersion 
and clumping feature is’ reversible as shown. by, 
It was Interesting to observe that 
this pH series of figure 7 maintained its 
features over a period of several months. The 
marked improvement in chemical durability of this 
glass in the H.SQ, solutions of concentrations vreatel 
than 50 percent (fig. 3 noted for other 
glasses, particularly in the repression of swelling of 
Corning O15 eleetrode glass in high concentrations of 
electrolytes 11:3] 

It will be of interest to learn whether fibers pro- 
duced from more durable silicate glasses show this 
same critical pH range for dispersal 
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essential 
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3.6. Thermal Expansion Characteristics 


Figure 8 shows the expansion curves for the glass 
from which the fibers are made compared with th 
expansion characteristics of Pyrex 7740 and Kimbl 
N 51-A. This figure gives an indication of the uppel 
temperature at which paper from this fiber 
glass can be successfully used, i. e., approximatels 
700° C for prolonged periods of time and somewhat 
higher temperature for shorter eveles Its deform- 
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ing temperature as determined by the interferometer 
is appreciably higher than either of the other two 
glasses shown. Inserted in the figure is a portion of 
the expansion curve for the fiber glass in an unan- 
nealed quenched condition. It shows the charac- 
teristic dip in the expansion as the temperature 


approaches the “critical range’ [20, p. 17] pear 
630° C, 
Early experiments had shown that the tensile 


strength of an experimental batch of paper, in which 
bentonite had been incorporated with the glass fibers 
was increased from 386 to 820 psi by heating to 
600° C, and that the density was also appreciably 
increased [2]. An adverse effect observed was an 
increase in brittleness of the fibers, a feature char- 
acteristic of glass fibers that have been leached and 
heated. These observations can be partially cor- 
related with the temperatures indicated by the 
expansion curve, which shows that the temperature 
of 600° C is within the lower limits of the “annealing 
range’ for this glass. The additional strength may 
have been brought about by initial stages of fluxing 
the fibers to each other or to the bentonite filler. 
For other glasses, heating after leaching is known to 
improve the chemical durability, decrease the hygro- 
scopicity of the surface layer, and greatly lower the 
voltage response of glass electrodes to hydrogen ion 
activity of aqueous solutions [21] The increase in 
density of the paper at 600° C would be accompanied 
by normal increase in density that is observed to take 
place in upon being altered from the 


“quenched” to the annealed condition [22] 


classes 


4. Relation of Observed Properties to Paper 
Quality 


The weakness of the fiber glass paper fabricated at 
values lower than pH 2.9 coincides with a sharp 
decline in the chemical durability of the glass The 
decrease in strength of the paper when formed at 
pH values higher than pH 2.9 is associated with a 
conspicuous tendency for the fibers to clump. This 
clumping probably interferes with the most advan- 
tageous orientation of the fibers for proper matting, 
and attainment of optimum strength of the finished 
paper, The evidence furnished by heterogeneous 
equilibria data and interference colors on the sur- 
faces of the interferometer test specimens, seems to 
support the assumption that the gelatinous silica-rich 
laver contributes to the resulting strength of the 
paper by causing the fibers to adhere at their points 
of contact when the paper is produced at the optimum 
acidity near pH 2.9. This feature, however, could 
be studied to advantage with other glasses of dif- 
ferent compositions. As to the ultimate tensile 
strength that might be attained for papers from glass 
fibers, it must be borne in mind that such values as 
400 psi are but a small fraction of the initial strength 
of the individual fibers. At the time of their produc- 
tion the tensile strength of glass fibers has often 
exceeded 2,000,000 psi [19, p. 72]. 

That the region near p!! 3 for maximum strength 
and quality of the paper from glass fiber is a range of 


pronounced transition of properties for typical silicate 


glasses can be readily visualized by inspection of 
figure 9. In this figure typical tensile strength data 
for paper from glass fiber are superimposed over the 
chemical durability curves for three silicate glasses 
(fiber glass, BaC 572, and Corning 015), typical 
voltage departure curve for the glass electrode, and 
a basic silicate-acid titration curve. All of these 
show a gentle to sharp trapsition in the upproximate 
region of pH 3. The same pH region is further 
emphasized by the peptization-pH_ series, (fig. 7 
and the position of the appearance of interference 
colors from reflection reducing films on the surface 
of the chemical durability test specimens, (fig. 3). 

There are other observations from this investiga- 
tion that are pot primarily concerned with the! 
production of paper from glass fibers, but that are of 
more than passing interest The data for the three, 
optical glasses given in table 2 and figure 1 when 
compared with results published vears earlier [4] 
show the peed for improvement in the chemical 
durability of these glasses; their poor durabilities have 
been contributory to the excessive deterioration of 
optical instruments and spectacle lenses under con- 
ditions of high atmospheric temperature and relative 
humidity [23]. Of equal interest is the possibility 
that peptization of glass fibers and paper strength as 
a function of pH may constitute useful tools in the 
study of the so-called “isoelectric pomt’” reported for 
silicate solutions near pH 2.9 [24, 25] 
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Adsorption of Nitrogen on Carbon Adsorbents at Low 
Pressures between 69° and 90° K 


Juan de Dios Lopez-Gonzalez,' Frank G. Carpenter,' and Victor R. Deitz 


Nitrogen adsorption isotherms were determined on three carbon adsorbents (two bone 
chars and an activated coconut shell charcoal) at three temperatures between 69° and 90° K 
and at relative pressures from 10~' to 0.9 \ ervostat was used to obtain temperatures 
below the normal boiling point of liquid nitrogen. The attainment of steady pressure 


readings at very low pressures required considerably more time than at higher pressures 
Corrections for thermomolecular diffusion were found to be appreciable at pressures below 
one millimeter of mercury 

At relative pressures below 10~4, the amount of nitrogen adsorbed depended upon the 
condition of the surface in regard to previous adsorption of nitrogen. A technique involving 
a series of Consecutive adsorptions at low temperatures and outgassings at room temperature 
led to a reproducible isotherm at low pressures. These data could be represented by a 
Freundlich equation. The Brunauer, Emmett, Teller surface areas of the carbon adsorbents 
were also determined 


1. Introduction 2. Experimental Procedure 
During the last deeade the main emphasis of A four-stage McLeod gage,’ shown in figure 1, was 


research in gas adsorption has been on the evalua- | added to a conventional adsorption apparatus. The 
tion of surface area. The success attained has been | calibration was made gravimetrically with mercury 
due largely to the development of the Brunauer, | and by pressure-volume data, using helium. Pres- 
Emmett, Teller (B.E.T.) theory of multimolecular | sures down to about 0.1 mm Hg could be measured 
adsorption and to the constructive critical work | with an accuracy of about 1 percent. An accuracy of 


, The fae : *.T. equa- Ea aEe 
that followed The fact that the BE. equa ‘Acknowledgment is made to W. V. Loebenstein of the Surface Chemistry 








Hon Was inapplicable to data at relative pressures Section for the design and fabrication of this gage 
below 0.05 was of little consequence because adsorp- 

tion measurements in this region were not needed {HIGH VACUUM 
in the surface-area determinations. Experimental i 


difficulties have been an additional cause for the 
scarcity of data on physical adsorption at low rela- 
tive pressures 

Adsorption data at very low pressures are neces- 
sary for the calculation of thermodynamic functions 
of the adsorbed phase. These functions are of spe- 
cial interest in the range of \ Cry low surface-coverage 
Ij. The interest in such adsorption isotherms has 
been further stimulated by the possibility of obtain- 
ing detailed information concerning the heterogeneity 
of a surface. The recent papers of Hill [2], Halsey 
3], Rhodin [4], Halsey and Taylor [5], Sips [6], and 
Cook, Pack, and Oblad {7] indicate some promising 
approaches. It is believed that surface heterogene1- 
ties and impurities are the more common cause of 
the variable properties of commercial carbon ad- 


sorbents. 
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Glevsteen [8S] for a number of carbon adsorbents. — REFRIGERANT 
; oR 
As these samples were still available, the opportunity CRYOSTAT 
was taken to extend the data to relative pressures 
of about 10°°. This corresponds to absolute pres- 
sures of about one-tenth micron of Hg. 
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2 percent was attained for pressures down to 107° mm 
Hg, 3 percent for pressures to 10-* mm, and 5 per- 
cent or more for pressures as low as 10°! mm, the 
lowest pressure attempted. The accuracy of the 


high-pressure (above 10 mm) readings was about 
0.05 mm Hg. 

Nitrogen from a evlinder of compressed gas (water- 
pumped was purified by passage over copper al 3500 
(* for the removal of residual oxygen and then 
through Ascarite and phosphorous pentoxide on glass 
to remove carbon dioxide and water vapor 
temperature of the 
sorption chamber was determined 
nitrogen or oxygen Vapor-pressure 
The temperatures 90.0 K were 
with liquid oxvgen and liquid nitrogen, respectively, 
boiling at atmospheric pressure. The temperatures 
72.2° and 69.5° K were obtained with a= special 
ervostat described in the following paragraph. In 
all cases the temperature of the adsorption chamber 


was constant to within about 0.05° K 


cotton 

The bath surrounding the ad- 
with either a 
thermometer 


and 77.7 obtained 


Cryostat 
The ervostat. see figure 2. consisted of a wide- 
mouthed 1-liter Dewar fitted with a rubber stopper 
and a 3.8-em inside diameter brass tube closed at 
one end, that reached nearly to the bottom. Liquid 


nitrogen placed in the annular between the 
brass tube and the Dewar was boiled under a fixed 
reduced pressure to obtain the desired temperature. 
In operation, air condensed inside the brass tube to 
form an excellent heat-transfer medium between the 
adsorption tube and tube. The vaporized 
nitrogen was drawn out of the ervostat through two 
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and contro 


The pressure controller was similar to that used ip 
vacuum distillations and consisted of a short mere ury 
manometer with electrodes, as shown in figure 9 
To start operations the stopcock on the controller 
was opened to the vacuum. When the desired] 
operating pressure was reached, the stopcock was 
closed and the controller thereafter maintained this 
pressure, As a considerable volume of nitrogen was 
vaporized, the piping, solenoid valve, and the vacuum 
pump were adjusted to provide an adequate capacity 
In oper: ation, the solenoid valve eveled about once 
every 2 to 5 sees Hach opening of the valve brought 
a sudden pressure but with a sufficienely 
large surge tank these flue tus itions resulted in varia- 
tions of only about 5-mm Hg pressure over the liquid 


o-oo 


de ‘crease 


nitrogen in the ervostat. The slope of the Vapor- ? 
pressure curve mn the operating region Is about 50} 
mm /deg, hence, the pressure variation corresponded , 


These fluetua- 
capacity of the| 


to changes of only about 0.1 deg C 
tions were smoothed by the beat 
apparatus, and variations of only +0.05 deg K were 
indicated by the vapor-pressure thermometer 

As a constant temperature was maintained over a, 
number of hours, it was necessary to replenish the 
supply of liquid nitrogen in the eryostat. A valve’ 
designed to operate when immersed in liquid nitro- 
gen is shown in figure 2. It consisted of a Pyrex ball 
joint without with the tubing to the ball 
joint sealed off and extended out of the liquid, where 
it served as a handle. In the off-position the ball 
was almost disengaged from the socket. In this 
position the opening in the socket was effectively 
closed by the shoulder of the ball. To open the 
valve, the bandle was twisted slightly to further 
disengage the joint. Despite minor and in- 
frequent clogging with ice crystals, the operation was 
quite satisfactory. 

The entire cryostat was lagged except for a small 
window to observe the level of liquid nitrogen 
About 3 liters of liquid nitrogen was consumed per 
hour in maintaining a constant temperature of about 
70° WK 


YTease, 


leaks 


2.2. Corrections for Thermomolecular Diffusion 


As there Is an appreciable temperature difference 
between the manometer and the adsorbent, there is 
the possibility of introducing appreciable error from. 
thermomolecular diffusion when measurements are} 


= 


made at very low absolute pressures A careful 
review - the available literature was made [9,10,11 
12], and the following equation of Liang [11] was} 


considered the most appropriate for calculating the, 


thermal transpiration factors 


. . ; 
R P On, p,\ . Du. p,.\ R,, 
Ons ~,.\)° | Du. PY | 
where 7’, is the observed pressure, 7’, the corrected 
pressure, and /? the correction factor NV is the prod-} 


uct of DP, by d, where d is the diameter of the tube 
connecting the cold to the warm part of the appa | 
ratus. FP 7, 7,)''", where 7, = temperature of thef 
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cold end of the connective tube (adsorbent tempera- 
ture) and 7; temperature of the warm end of the 
connective tube (room temperature). The terms 
ye ANd Bye are empirical constants: (1) aye=2.52 
independent of temperature and (2) Bye=7.68 
;—f,,), the pressure being measured in millimeters 
of mercury and the tube diameter in millimeters. 
The factor ®, is the empirical relative ‘“pressure- 
shifting factor’? and is dependent upon the gas but 
not the temperature, It is related to the collision 
diameter the molecules. With the above 
listed values of @ and 8, ®y.=—1 and @y =3.28. 


of cas 

lt is evident that as the pressure approaches zero 
(Y-0); then PR approaches its theoretical minimum 
value, ?,,. Also, for high pressures (large values of 
X), R approaches unity. In the present work, the 
connecting tube was 1.0-mm diam, and nitrogen was 
the gas used at low pressures. In figure 3 are shown 
the variations of 7/2 with pressure for the present 
ease in the interval from 0.0001 to 1 mm Hg for the 
four different temperatures. The values of F for 
pressures above | mm Hg are very close to unity. 
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2.3. Materials 


The three carbon adsorbents were among those 
previously investigated [8]. Char 2 was a new bone 
char obtained from the manufacturer. Char 27 had 
been employed in a cane-sugar refinery for a con- 
siderable time and had been removed from the char 
stock because of its high specific density (1.3 g/ml 
compared to 0.63 g/ml of char 2). The sample of 
char 27 used in this work differed from that previ- 
ously reported [8] in that foreign substances visually 
detectable were mechanically removed. Char 1, an 
activated coconut-shell charcoal, was identical with 
that previously reported [8]. 


3. Experimental Procedure 


The three adsorbents were simultaneously out- 
gassed at 400° C for 16 hr at the beginning of the 
experiments. The pressure over the freshly out- 
gassed adsorbent before any addition of nitrogen 
was always of the order of 10-° mm Hg. In all 
cases, the first addition of nitrogen was very small, 
so that equilibrium was attained at a very low 
pressure. Then, more nitrogen was added, and the 
volumes adsorbed were obtained at increasing equi- 
librium pressures. After several points had been 
determined, the samples were outgassed for 16 hr at 
room temperature. This process of adsorption and 
outgassing at room temperature was repeated Ssev- 
eral times for each sample. Because of the great 
experimental difficulty in obtaining desorption data 
in the low pressure range, all of the results reported 
here are for adsorption. 


4. Results 


4.1. Attainment of Steady-State Pressure 


The time required to reach a steady state was 
found to be a function of the final pressure of the 
adsorbed nitrogen, and depended also upon the stage 
of outgassing operations. After the first outgassing 
at 400° C, the steady state at the lowest pressure 
reading was attained in 1 to 2 hr, but in the meas- 
urements after the final outgassings at room tem- 
perature the attainment of a steady state required 
a much longer time, as is illustrated in figure 4. In 
all cases, the attainment of a steady state after the 
introduction of the first portion of nitrogen required 
the longest time. Curve A in figure 4 was observed 
for the lowest point of the isotherm for char 27 at 
90° K (see fig. 6) and a steady-state pressure of 0.4 u 
was realized after 10 hr. In several other similar 
cases it was observed that such a steady state per- 
sisted after 13 additional hr. Curve B of figure 4 
illustrates that only 3 hr was required when the 
steady-state pressure was about 200 uw. When the 
pressure was in excess of 1 mm, the steady state 
was reached in about 20 min, in agreement with 
earlier experience [13] (curves C, D, E, and F, fig. 4). 

The influence of consecutive outgassings on the 
quantity of nitrogen adsorbed was of particular in- 


terest. Figure 5 contains the low-pressure adsorp- 
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example that the amount adsorbed in- 


creased with consecutive outgassings. By interpo 


lating each curve to an equilibrium relative pressure 


of 2 10 


the volume of nitrogen (milliliter STP 


adsorbed was 0.004 ml after the initial outgassing at 
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at room temperature. The isotherms become coin- 
; 
cident above a relative pressure Of about 10-7, in 


which recion this effect was not evident 

In the first adsorption experiments with chars 2 
and 27 after the first outgassing at 400° (. a small 
Increase In pressure was observed in the svstem when{ 
the first portion of nitrogen at very low pressure was 


introduced into the adsorption chamber This in-! 


crease Was observed only during the first few minutes f 
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of contact, and then the pressure decreased in a nor- 
mal way toward equilibrium. When this effect was 
observed, the final values of the relative pressures 
were notably high. A possible explanation is pre- 
sented in a later paragraph. 

The times required to realize a steady state for 
char 27 were only about half those required for chars 
2 and 1 at the same final pressure. Periods of about 
50 hr were necessary to attain equilibrium in chars 
| and 2 at the lowest pressures obtained. 


The volumes of nitrogen adsorbed for char 2 
(69.5°. 77.7°, 90.0° K), char 27 (72.2°, 77.77, 90.0 
K). and ehar 1 (77.7°, 90.0° IX) are given in tables 
1, 2, and 3, respectively The data recorded are the 


final steady-state values observed after a number of 
consecutive temperatures suffi- 
cient to vield reproducible values for the volumes of 


outgassings at room 


nitrogen adsorbed 
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4.2. Calculation of Isotherm Equation 


An isotherm equation was sought, valid for the 
region of the lowest measurable pressures, which 
could then be used with some confidence to extra- 
polate to the volume adsorbed at pressures too low to 
measure. Such an extrapolation is necessary in the 
calculation of thermodynamic functions of adsorbed 
molecules [14]. 

Plots of the volume adsorbed as a function of 
relative pressure for chars 2, 27, and 1 vielded curves 
that became steeper as the pressure was reduced, as 
shown in figure 6 for char 27. A linear behavior 
was found, however, on a Freundlich plot for all 
temperatures and all chars from the lowest relative 
pressures up to about 107%. Typical plots are 
shown in figure 7. The isotherms observed at the 
lower temperatures were linear up to a higher value 
of relative pressure. For example, in char 27 at 


72.2° K a straight line was observed up to p/p,.™3> 
10-*; at 77.7° K the straight line was valid up to 
5107, and at 90.0° K up to 4X10~*. A similar 


was observed in chars 1 and 2. Parameters 
Freundlich equation were determined, and 
given in table 4. 


behavior 
n of the 
these are 


TARLE 4 lsotherm constants for adsor ption of nilroge nm 
B.E.T. parameters 
I I Freundlich 
{ 
ire constant 
| c Fa Surtace area 
A nl m?/¢ 
3. 92 397 4 4. 385 1,744 
wt s2 47s t 1. 601 1.741 
{ Hu 7.40 28. 65 104 4.2 121.9 
, 77.9 7.39 97. 33 I8v 4, 383 119.8 
| Ww l 25. 54 183 4. 601 117.5 
2.2 7. 30 6. 95 110 4. 205 20. 85 
7 7. 30 6. 60 123 4. 383 93 
00.0 7. 30 6.28 107 4.601 2s RU 
Pack based on density of liquid nitrogen for converting V, to 
rface area 
Calculated from the Langmuir isotherm 
* 
| 
} 
' 
=x 8 x x 6x 5 


FiGuRE 6 Vitrogen adsorption isotherms on char 27 at 72.2 


~ 6°. and 90° K. 
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4.3. Surface Area Calculations 


The adsorption data for chars 2 and 27 were found 
to give linear B. E. T 
pressures between 0.05 and 0.35 
and ¢ were evaluated from the B 
a free surface 


The values of V 


ri 
4 
4 


p l ( | 
V(p Pp) eV eV, P BP i 
s the volume adsorbed 
the saturation 


(reported at STP 
P, and i 


where V 
at pressure 


Is 


pressure 


The values of the constants V,, and ¢ are given in 
table 4. The corresponding surface areas were 
based on a packing factor calculated from the 


density of liquid nitrogen [8]. The values given in 
table 4 were calculated from the adsorption data 
recorded in tables 1, 2, and 3. V’,,, also calculated 
from the data obtained in the first adsorption run 
after outgassing at 400° C, was found 26.5 
compared to 25.5 for the steady-state value. De- 
spite the large difference between the first portion of 
this isotherm and the corresponding part of the 
isotherm obtained after several outgassings at room 


to be 


temperature, there was close agreement in the values 


of V,, obtained from the two sets of data 

The adsorption data for char 1 did not give a 
linear B. E. T. plot. The correct values of V,, for 
char | are given by Langmuir plots, which were 


found valid in the relative pressure range between 
0.2 and 0.8. The values of V,, at 77.7° and 90.0° K 
are also given in table 4 


5. Discussion 


5.1. Attainment of Steady-State Pressures 


As the carbon adsorbents investigated were essen- 
tially pyrolysis residues of shell 
and not subject to a defined state of decomposition, 
it is highly probable that some surface complex 
containing carbon, hydrogen, and oxygen remained 
with the This could decompose into 


bone or coconut 


adsorbent. 


plots in the range of relative 


equation for 
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ee 
f the nitrogen adsorptior 
at 72.2 6°. and 90° K 
CO,, CO, and H, by continued thermal treatment 


Some of this complex could be removed during the 
outgassing at 400° C, but an appreciable portion 
always remained in the adsorbent. The 
of this residual complex could have a major influence 
on the observed steady-state pressure in the range 
of very low pressures 


presence 


The experimental results suggest that one of the 
factors responsible for the behavior observed at low 
pressures may be directly related the residual 
complex. The first portion of nitrogen introduced 
was able to displace some constituent of the complex, 
and, asa result, the liberated decomposition products 
such as ( | ). (( ),, and I] . contributed to the observed 
pressure Unfortunately, means were not availabl 
to analyze the products of outgassing at the low 
pressures emploved, but this would be Pa | desirable 
future experiment. The consecutive outgassings at 
room temperature, together with introduction 
more nitrogen, served to remove the complex until 
a stage of relative stability was reached, as shown 
by the final reproducible values. The particular 
program of outgassing is obviously one of a possibl 
large number. Any alteration in the technique of 
outgassing, such as a heating to 400° C, 
would no doubt have had an additional influence on 
the number of adsorption sites shared between the 
surface complex and the adsorbed nitrogen 


to 


of 


second 


Another factor of influence is associated with th 
diffusion of gas molecules at low pressures. Becaus¢ 
of the intricate structure of the porous material, more 
time is required for the nitrogen to locate the adsorp- 
tion sites of minimum potential energy. A similar 
argument may be valid in connection with surface- 
diffusion phenomena. Some such consideration is 
required to explain the fact that even relatively well 
outgassed adsorbents (400° C for 16 hr required long 
periods to attain equilibrium in the low-pressure 
range {15}. 

Other things being equal, relatively shorter periods 
> hed 


were required to attain a steady state for char 
compared to char 2. It has been shown elsewhere 


[16] that the crystalline structure of basic calcium | 


phosphate (about 90 percent in these materials) was 
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Was 


more highly developed in char 27. In this case the 
contribution of the voids in the porous structure to 
the total adsorbing surface was less. In addition, the 
carbonaceous residue in char 27 had a greater stability 
as a result of many heatings to which it had been 
subjected in the commercial revivification process. 


5.2. Validity of Isotherm at Very Low Pressures 


The experimental data show that a Freundlich 
isotherm is valid in the measurable range at very low 
pressures. When a Langmuir isotherm is valid at 
low pressures, Henry’s law, V=az, is automatically 
obtained as the limiting case. Hill [17] has shown 
from statistical mechanical-thermodynamic —con- 
siderations that all adsorption isotherms must ap- 
proach Henry’s law. Hence, in the present case 
where Freundlich isotherms are observed at the 
lowest pressures, it is necessary to postulate a transi- 
tion range below which one must realize the linear 
dependence of adsorbed volume on pressure. It 
not known where this range exists for the ad- 
sorbents investigated. 

For Henry ’s law to be obeved, the slope of the 
plots in figure 7 must be unity. It may be seen from 
these graphs that the first observed points at the 
lowest pressure for the isotherms at 90.0° and 77.7° K 
are displaced in the correct direction. However, it 
will be necessary to obtain additional data at even 
lower pressure to be sure of this behavior. 
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5.3. Surface Area 


When the ratio V/V, was plotted (see fig. 8) as a 
function of p/po from 1077 to 10~‘, ii was found that 
the curves for chars 2 and 27 were in coincidence at 
90° K. In this region no more than 40 percent of the 
surface was covered. This agreed with previous 
work in the range 0.1 to 0.3 [8], where the entire 
surface was covered. However, at the lower tem- 
peratures, the curves for different adsorbents were 
separated to a significant extent, which implies a 
different interaction of the nitrogen molecules with 
the surface. For the case of a given adsorbent at 
2 or 3 temperatures, the curves were separated to 
greater extents than indicated by the corresponding 
changes in the densities of liquid nitrogen. The 
values of V/V, are largest at the lowest tempera- 
ture for a given value of p/ppo. 

The surface areas given in table 4 indicate an 
increase with decrease in the temperature. It is not 
known whether or not this is within the experimental 
error of the VV, determination. It may also be 
related either with unknown changes in density of 
the adsorbed nitrogen phase or with the possibility 
that the effective van der Waal’s diameter of the 
adsorbed nitrogen molecule at the different tempera- 
tures is critical with the dimensions of the voids. 
Therefore, the surface available to a nitrogen mole- 
cule, especially in porous materials, might be corre- 
spondingly greater. 


5.4. Conclusions 


The attainment of equilibrium at very low pres- 
sures is a complicated process, in which the nature of 
the surface and the impurities associated with the 
solid surface play an important role. Variations in 
the outgassing techniques can change this part of 
the isotherm to a notable extent. 

As the nature of the predominant forces respon- 
sible for adsorption changes significantly with surface 
coverage, no known isotherm equation can be ex- 
pected to be valid from r=0tor=1. The Freundlich 
equation appears to be valid at very low pressures in 
a great many cases, but it is not valid when the 
surface coverage approaches a monolayer where z is 
approximately 0.1. The B.E.T. equation appears to 
be frequently valid between the point where the 
Freundlich relationship fails and relative pressures of 
about 0.35. In this region (0.05 to 0.35) the factors 
that are not taken into account in the B.E.T. model 
would not contribute appreciably to the volume ad- 
sorbed. The behavior above 0.35 in relative pressure 
has at present no adequate theoretical basis. In 
general, one can say that the nature of the pre- 
dominant forces involved in any adsorption process 
is a function not only of the nature of the adsorbate, 
adsorbent, and temperature, but also of the pressure 
range, which determines the fraction of the surface 


covered, 
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| Stress-Strain Relationships in Yarns Subjected to 
Rapid Impact Loading: 3. Effect of Wave 


{ Propagation 
( Jack C. Smith, Frank L. McCrackin, and Herbert F. Schiefer 
The tensile behavior of a Hookean material, elongated by rapid impact at one end 
has been calculated, using a theory in which wave propagation is considered. As a result 
of these calculations, limits have been established on the applicability of a simpler theory, 
’ discussed in a preceding publication, in which wave propagation was neglected. 
1. Introduction 
When a filament of uniform cross section is clamped at each end and rapidly elongated, 
i data for a stress-strain curve can be obtained by measuring the force at one end and the elonga- 
: tion of the whole filament, both as functions of the time. In this case the assumption is made 
that stresses and strains in the filament are uniform along the length and can be represented 
; by the stresses measured at an end and by the average elongations. This assumption is no 
longer valid when the rate at which the filament is lengthened approaches in order of magnitude 
the velocity of propagation of a tension wave along the filament. 
In previous work * this effect of wave propagation was disregarded. It is the purpose 
here to determine how previously obtained results must be reinterpreted as a result of the 
wave-propagation effect. 
2. Wave Theory 

Consider a filament ving along the negative z-axis, with one end at the origin and having 
the other end fixed at the point a lL. At time f=0 and subsequently, the end at r=0 is 
constrained to move with velocity 7 in the positive direction. Assume that the tensile behavior 
is in accordance with Hooke’s law. Then, after impact, a strain pulse of magnitude €)==7/e 
is propagated down the filament with a velocity ¢. When the pulse is reflected at the fixed 

; end, the strain is increased to 26. A similar increase takes place at each reflection. 

This behavior is shown in figure 1, in which the strains at several points along the filament 
are plotted as functions of the time. The strains are thus seen to vary discontinuously with 
time. These discontinuous increases in the strain become more pronounced as v is Increased. 

If we require that « be less than 1 percent, in order that wave-propagation effects in most 


textile fibers may be safely disregarded, we must limit ¢) to values less than 0.01 ¢, which in the 


case of high-tenacity nylon is approximately 25 m/sec 
In the present experimental method a filament is terminated at s L by a mass nw, 
where » is a number and w is the mass of the filament. At x0 the filament is terminated 
by a mass so large that its momentum will be practically unchanged by reaction forces from the 
filament and attached mass. At time ¢=0 the head mass is given a velocity ¢. The behavior 
of the system is recorded photographically... The data of interest are the elongation of the 
filament and the reaction force applied to the filament at the tail mass, both as functions of 
time. From these data the experimental stress-strain curve is constructed. 
‘ kw I } the Office the Quartermaster General, Department of the Army ind was presented in part at the November 
‘ es RK t 
Frank L. MeCrack I hert F, Schiefer, Jack C. Smith, and Walter K. Stone, Stress-strain relationships in yarns subjected to rapid 
ling: 2. Breakin elocit ie train energies, and theory neglecting wave propagation, J. Research NBS 54, 277 (1955) RP2590 
W rkK.s I F. Sch ind George Fox, Stress-strain relationships in yarns subjected to rapid impact loading: 1. Equipment, 
tua 1d h NBS 54, 269 (1955) RP2589 
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The theoretical behavior of this system, W hen Hooke's law applies,’ Is most easily calculated 
by finding suitable solutions of the wave equation 
Oru Oru 
of Oy 
In this equation, ¢ is the time variable, and « denotes the position of a cross section of the 
filament relative to fixed coordinates, when the filament is in the unstrained state; u(z,f) 
denotes the distance the cross section has moved from its original position at sz, and ¢ is the 
velocity of propagation of an elastic-tension wave along the filament. This velocity is given 
by c= Ep, where F is Young’s modulus of elasticity for the unstrained filament, and p is 
the density of the unstrained filament material 
These solutions must satisfy the boundary condition at the head: 


u(Q, tf) Pot (2 


A suitable solution (for the unreflected wave) is given by 


(t+") ( t+" 0) 


e(r,t 0 ( t+" <0) 


where u represents a wave of velocity ¢ traveling along the filament in the minus z-direction, 


filament 


at Ir () 


and ¢ is the local strain, or increase in length per unit length, set up in the filament in the wake 


of the wave It should be noted that both ~ and ¢ are equal to zero for values of ¢ ric, ie 


no strain or displacement is present at a point in advance of the wave. 


When the wave arrives at the end of the filament it is reflected and travels back in the 


positive direction. The strain between the reflected wave front and the tail mass is represented 





‘ The stress-strain behavior of most textile fiber ribed by Hooke iw or t all strain At larger strait the solutior 
lerived here provide only rough approximations to the behavior of an actual textile fiber but are adequate for the purp ! ry 

In order to take into account the significant redu n irea of the he f ‘il ‘ | ! Sst | 
that Hooke’s law is obeyed in such a way that the product of Young’s modulus and ere ectional area remains constant Che small effect of n 

tion associated with the lateral contraction ha 
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> Figure 2. Strain-time curves for a filament loaded at FiGuRE 3. Stress-strain curves for a filament loaded 
the tail with a mass 10w and rapidly elongated at with tail mass 10w and rapidly elongated at the head 


the head at velocity at velocity vo 


as the sum of the strains caused by the incident and reflected waves. If the tail end of the 
filament were fixed, this strain would be 26. However, as the tail mass can move, this strain, 
initially 2e at reflection, will decrease with time. The rate of decrease depends upon the tail 
mass and upon the order of the reflection; i. e., whether it is the first or a subsequent reflection. 
Reflections also occur at the head mass. <A different solution of the wave equation is required 
to describe the state of the filament after each reflection. The first 10 of these solutions are 
tabulated in appendix 1 

By using these solutions it is possible to calculate curves of average strain versus time, and 
local strain versus time. These curves are given in figure 2 for the case where the tail has a 
mass of 10w (n= 10 The curve of local strain, e(—L,?¢), versus time is jagged, the decrease 
in « between successive reflections becoming greater with each reflection. 

As the force at the tail, which is found experimentally from the acceleration of the tail mass, 
is proportional to the strain at the tail, it is possible to construct a theoretical stress-strain 
curve of local tail stress (strain) versus average strain. Such a curve is given in figure 3. 
Experimental curves are similar to this curve in that they show dips in the vicinity of maximum 
stress. There is also a delay in the initial rise of stress by the time it takes for the tension wave 


; to be propagated to the tail mass 


3. Theory, Neglecting Wave Propagation 


A solution, applicable at low testing speeds, can be obtained by disregarding wave propaga- 
tion along the filament, and assuming the mass of the filament to be concentrated partly at the 


tail and partly at the head. The resulting solution, derived in appendix 2, is 


A 
uel | @d— yn’ sin (4) 
\ i] 
€ = . fe 
yr sin ; (0) 
€q yi 
where 6=— (e/L)t. e=—rm/e. n’=n-+f, and fis the fraction of the filament mass assumed concen- 


trated at the tail 

It is helpful to think of @ as time expressed in L/e units. Likewise, one can think of [e/e] as 
strain expressed in « units. Here « and ¢ have only formal significance in order to put the 
solution in a form analogous to that obtained when wave propagation is considered. 
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4. Comparison of Wave Theory With Theory Neglecting Wave Propagation 
and 11 are shown in figure 4. Values of ; 


10, 1043, 
10. using the formula 


Plots of le € versus @ for values of n’ 
é/e, have also been computed from the wave theory for the ease n 


€ f | c 
A u | 6 
€ é L Lr ’ 
These points are plotted as circles and are seen to fit best to the curve drawn for n’= 104), 
This indicates that for a mass ratio of 10, the nonwave theory results agree most closely with > 


those for the wave theory when half the mass of the filament Is assumed concentrated at each 


end (f s). 
According to the nonwave theory, the maximum average strain that can be obtained in a 
; 


filament is given by 
. 
[2] ante , 
€ ’ , 


and this strain ts obtained when @ 


This suggests that the maximum average strain €/€ as computed by the wave theory ts approxi- 
mately equal to ,+4 and occurs for a value of 6, approximately equal to 2 ,n+3. That 
this is a good approximation is demonstrated in table | 
According to the nonwave theory (appendix 2 
> 
+ 


ale. 


oY 1 OL (li L log 
Vo. 


where v, is the lowest impact velocity to cause rupture in a filament to which a tail mass, 7” 
is attached, o@ is the tensile stress expressed as force divided by the cross-sectional area of the 


unstrained filament, and e€, is the strain at which rupture occurs 
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If the values of é, and o are independent of the rate of elongation, then log \ «| p— ale 
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When log (n- b) is zero, we have the extrapolated velocity 


\ ade. (10) 
P Jo 


When Hooke’s law applies, we have 


According to the wave theory, , //p is the velocity of propagation, c, of an elastic tension wave 
along the filament, and thus 


the velocity just sufficient to cause immediate rupture at 2=0. 
We call the velocity defined by (10) the “limiting breaking velocity”. The concept of a 
limiting breaking velocity has been discussed, in a previous paper (see footnote 2). 


As an example, consider a Hookean material for which ¢€,=0.12, and e=2,500 m/sec. 
From 7), for the cause of rupture, we have €, fey la , 0,/C)v~nr-7 ; Or 0. 


In figure 5 are plotted values of ¢, computed from these assumed values of e, andc. The 
result is a straight line giving an extrapolated value of 300 m/see as the limiting breaking 
velocits 

The breaking velocities have also been calculated on the basis of the more exact theory in 
which wave motion is considered and Hooke’s law is assumed. Average strains are calculated 
from (6), and the breaking velocities are then computed from the formula 


Ce, 


, m= ° (14) 
€(/i) 
e max 


These results are also plotted in figure 5. It is seen that these breaking velocities, for a 
Hookean material of constant e,, are lower (for small n’s) than those predicted by the nonwave 
theory. 

In actuality the strain at a given point along the filament does not change continuously 
with time but increases suddenly by an amount €) when one of the strain pulses passes through 
the point. Rupture then occurs whenever the local strain exceeds the rupture strain. 
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In figure 6 are plotted local strains at the head, tail, and midpoint for a filament loaded 
with a tail mass 10w (n= 10 Note here that maximum strain occurs at the head on the third 
reflection (0 6). This corresponds to a lowest breaking velocity of 68.6 m/sec 


Increases in the local strain alwavs occur first upon reflection at either the head or the tail 
Increases in the local strain at an intervening pomnt occur later, but these strains never exceed 
the head or tail strains in magnitude. For this reason, it is to be expected that a uniformly 
strong filament will break only at the ends 

If the maximum values of local strains [€/€o]max are substituted into eq 14, the lowest 
breaking velocities, v,, for head and tail breaks can be computed These breaking velocities, 
also plotted in figure 5, are seen to be the lowest of all 

The curves of local strains versus time for an actual material, having approximately the 
properties chosen for this example, would not be as jagged as those shown in figure 6 for an 
ideal Hookean material. Thus, the velocities found experimentally for this actual material 
could be expected to lie between the bounds established in ficure 5 by the curve of the velocities 
computed from average strains in the filament, and the curve of velocities computed from the 
extreme local strains at the head or tail of the filament. The breaking velocities computed by 
any of the methods considered are the same when vw is greater than 50. It is thus possible, in 
this example, to determine a limiting breaking velocity by extrapolation of experimental data 
taken with » greater than 50, provided ¢, remains constant and Hooke’s law applies. This 
criterion also applies, provided the energy density at rupture for the material remains nearly 
constant at the impact velocities used for determination of the limiting breaking velocity. 

In the example just considered, the breaking velocity for »—50 is 42 m/sec. A tensile 
impact at this velocity produces a strain pulse of ¢.—1.6 percent. This suggests a criterion 
applicable to any textile fiber, for which the energy density at rupture remains nearly constant 
at high strain rates. For such a fiber a limiting breaking velocity may be found by extrapola- 
tion of experimental data taken such that all 7, used are less than 0.02 ¢, where ¢ is the velocity 


of propagation of the tensile strain pulse along the sample for the material tested 


5. Summary and Conclusions 


When the stress-strain characteristics of a textile fiber are measured at longitudinal impact 
rates exceeding 10 m/see the effects of tension-wave propagation along the fiber must be con- 
sidered. These initial and reflected waves cause the local strains to vary with time in a step- 
wise fashion. The average strain in the filament, which is the quantity measured experi- 
mentally, does not adequately represent the state of strain under these conditions. As the 
measured stress is roughly proportional to the local strain at an end, experimental stress-time 
curves also vary in a stepwise fashion. Stress-strain curves likewise exhibit fluctuations 

The tensile behavior of a Hookean material having a breaking elongation independent 
of testing speed has been calculated, using a theory in which wave propagation is considered 
The results of these calculations have been used to modify the predictions of a simpler theory 
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that neglects wave propagation. Thus, it has been shown that a filament weighted at the tail 
and pulled rapidly at the head at constant speed will attain a maximum average elongation 
proportional to ,n+4 after the tension wave has been reflected (a 2),n+4 times, where n is 
the ratio of the tail mass to the mass of the specimen. 

If logarithms of the lowest breaking velocities for each n are plotted versus log(n+ 4), the 
theory neglecting wave propagation predicts that a straight line of slope —4 is obtained, 
provided that the rupture-energy density is nearly constant at high testing speeds. An extra- 
polated velocity (for log(n+ 4) equals zero) is thus obtained, which is termed the “limiting 
breaking velocity”? of the material; i. e., the velocity above which the sample will always 
rupture at the head immediately upon impact. Both theories predict the same limiting break- 
ing velocity, provided extrapolation is made from data taken at impact speeds less than 2 
percent of the velocity of propagation of a tension-wave pulse along the material. — If this limita- 
tion is observed, the theory neglecting wave propagation is valid, and extrapolation to the 


limiting breaking velocity may be made. 
6. Appendix 1. Wave-Theory Solutions 


In table 2, vw gives the displacement due to the original incident wave, proceeding toward 
the tail mass. After reflection at the tail, the displacement, uv, is computed from the sum of 
vw», and u;, the displacement of the reflected wave. Similarly, after the first head reflection, 
displacement is computed from the sum of uw, %4,, and uw. Displacement due to each reflected 
wave is a function of a 2 parameter, which in turn is a function of z and ?¢. Only positive values 
of these 2 parameters are allowed. Negative values represent position-time events in advance 
of the wave, where u=9. 

Velocity functions obtained by differentiating the displacement functions with respect to 
time are tabulated in table 3. Local-strain functions, obtained by differentiating the displace- 
ment functions with respect to the position coordinate along the filament, are given in table 4. 

The method of obtaining these solutions is demonstrated in the following derivation for u,. 


The initial conditions that must be satisfied for the first reflection at the tail are at 


J te *$ > 2Z=0, 2,=0 


Displacement at the tail just before reflection equals the displacement at the tail just after 
reflection, or w(0)=up(0)+u,(0). Hence u,(0)=0. Velocity at the tail just before reflection 


equals the velocity at the tail just after reflection, or 0=7,(0)- r(0). Hence v, (0) Uo. In 
addition, we have the boundary condition at x 
O° ' O 
nw —, [uUo( 20) + 1(2)|2--.=EA [ o( Zo) + Uy (2))]-- —2 
ot Or 
ol 
re +. . Oo 
nw —; u,(z)=EKA-—+EA u, (21) 
of- c Or anal 
or 
ih ..t EKAd 
nw; u,=—EA uy 
dz ¢ Cc dz 


where 2—t— L/e. and A is the cross-sectional area of the unstrained filament. If we let 
nwe npLe L 
T ; ; yi ’ 


we obtain 
dru, l du, i 
dz T dz T 
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The solution of this equation that satisfies the initial conditions is 


If we replace (t—L/c) by Pn DEE Se c, we obtain the required solution for u 
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Appendix 2. Solution, Neglecting Wave Propagation 


Consider a filament of length 1, cross section A, and density p and having a mass n’pAL 
attached at the tail, where n’ is equal to n+-f, n is the ratio of the tail mass to the mass of the 
filament, and / is the fraction of the mass of the filament assumed concentrated at the tail. It 
is shown in section 4 that the best value for fis }. 

At zero time and subsequently, let the head end of the filament be constrained to move with 


velocity %. The equation of motion is 


P Pu ’ 
i pAL it F 15) 


where uw is the displacement of the tail mass, and Fis the accelerating force transmitted through 


the filament to the tail mass. 
The strain in the filament, ¢, is given by 


Nol u , 
a ” 
When Hooke’s law applies, 
F=AE¢, 17) 
where / is Young’s modulus for the filament material. Substituting (16) and (17) into (15) 
we obtain 
du 
det Wo = worl, 18) 
where 
ae 19) 
Lyn’ 
and 
E 
()) 


A solution of (18) satisfving the boundary conditions, u=0 and du /dt=0 at t=0, is 


| 
ut [ sin wt | 21) 


If we let ¢ Lie 6, and e¢ rnc, We obtain 
ers 4 
uel | @— yn’ sin 22) 
yf 
, 6 
€ €ésy fe sin )*3) 
y fi 


In a previous paper (see footnote 2) the following equation was derived: 


ii- 4 ei ! 


> >) ») oo 
*'+log ale, 24) 
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log v, } log 


is the lowest impact velocity sufficient to cause rupture in a filament to which a tail 


where ¢, 
mass 7 times that of the filament and a head mass m times that of the filament, are attached. 
For comparison with the work discussed here, we let m—> ©, obtaining 
log r, log (n- +) ' log _ ale 25) 
Vo. 
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Thermal Expansion and Phase Transformations 
of Low-Expanding Cobalt-Iron-Chromium Alloys 


Peter Hidnert and Richard K. Kirby 


Coefficients of linear thermal expansion of some cobalt-iron-chromium alloys are reported 


for various temperature ranges between 


6: 


and 
perature, chemical composition, heat treatment, etc., 


the effects due to tem- 
Some of the alloys 


t-S00° C, and 
were determined 


investigated have coefficients of expansion less than those for fused quartz and ordinary 


invar between 65° and + 60° C. 
have ¥ 
to y on heating to high temperatures 
and Ac, temperatures at about 600° C 


to 1,000° C 


Some of the low-expanding cobalt-iron-chromium alloys 
>a transformations on cooling to low temperatures, and the resulting a-phase reverted 
Ar, temperatures were observed as high as 

The effects of various heat treatments from 
on the transformations were investigated, and the resulting changes of thermal 


10° C 
196 


expansion were correlated with the structure of these alloys 


1. Introduction 


In 1934, Masumoto [1]' reported the results of his 
investigation on the linear thermal expansion of 
cobalt-iron-chromium alloys containing more than 
19 percent of cobalt and 5 to 20 percent of chromium. 
These alloys were prepared by mixing suitable pro- 
portions of electrolytic cobalt, electrolytic iron, and 
chromium (with 2.7% of impurities) and melting 
these metals in alumina crucibles in a hydrogen 
atmosphere. The melts were cast in iron molds 26 
em long and 0.6 ¢m in diameter. Samples 10 ¢m in 
length were cut from the ingots and heated for 1 hr 
at 1,000° C in a vacuum furnace and then cooled 
slowly in the furnace. The chemical compositions 
indicated by Masumoto for these alloys were not 
obtained by chemical analyses [2] but were computed 
by him from the initial charges. 

Some of these alloys investigated by Masumoto 
have very small coefficients of expansion. For 
example, the coefficient of expansion of an alloy 
with a nominal composition of iron 36.5 percent and 
chromium 9.5 percent is 0.1 10-°/deg C between 
20° and 60° CC. This expansion coefficient is less 
than those of such well-known low-expanding 
materials as fused quartz [3] and the usual grade of 
invar [4]. The investigation of these low-expanding 
cobalt-iron-chromimum alloys was made over the 
range from liquid-air temperature to about + 250° C 
Masumoto stated that the observations above room 
temperature were taken during beating, as the 
expansion was almost reverisble for heating and 
cooling, and the observations below room = temper- 
ature were taken during cooling. 

In his paper Masumoto also reported results of 
investigations some magnetic, electrical, and 
corrosion-resisting * properties of one of these low- 
expanding alloys. Data on Young’s modulus and 
the temperature coefficient of Young’s modulus for 
various cobalt-iron-chromium alloys were reported 
by Masumoto and Saito [5]. 


On 


Figures in brackets indicate the literature references at the end of this paper 
Chis alloy resisted corrosion in a dilute sodium chloride solution to a much 
extent than ind th Masumoto named it “Stainless-Invar.”’ 


eater nvar, erefore 
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During World War II, Eastman Kodak Co., 
Rochester, N. Y., obtained a Government contract 
from the Office of Scientific Research and Develop- 
ment to produce a low-expanding alloy having a 
coefficient of expansion less than that of a good 
grade of invar at atmospheric temperatures. After 
a search of the literature and some preliminary tests 
with the cooperation of other companies, Eastman 
Kodak Co. authorized the Unexcelled Manufacturing 
Co., Cambridge, Mass., to prepare some low-expand- 
ing cobalt-iron-chromium alloys similar to those 
investigated by Masumoto. Peter Hidnert and 
Wilmer Souder of the National Bureau of Standards 
made determinations of linear thermal expansion on 
a number of these allovs over the range from ‘ 
to +60° CC. A restricted report of the results of this 
work was prepared by Edgar D. Seymour of Eastman 
Kodak Co. for the Office of Scientific Research and 
Development. The following summary was taken 
from an unclassified extract of the report: 


65 


From a search of the literature for a low-expanding 
metal, Masumoto’s cobalt, iron, and chromium alloy gave 
the most promise, if the melting, heat treatment, and 
analysis problems could be solved. We believe these 
problems have been solved and the technique has been 
perfected so that it is possible to specify the composition 
to 0.1% in chromium and iron, An alloy has been made 
with a coefficient of expansion of less than + 0.5 10-*/° C 
over the range from — 65° to + 60° C with the coefficient 
nearly zero from 20° to 45° C and also in the region of 

35° C. This is from two to three times better than 
invar obtainable on the market. 

The data reported by Seymour on thermal expan- 
sion of cobalt-iron-chromium alloys is included in the 
present paper. Additional data by the present 
authors on these alloys and other cobalt-iron- 
chromium alloys are given for various temperatures 
between —65° and +-950° C, 


2. Alloys Investigated 


The cobalt-iron-chromium alloys listed in table 1 
were prepared by John Wulff of the Unexcelled 
Manufacturing Co., Cambridge, Mass. The alloys 
were made from cobalt, iron, and chromium powders. 
A master alloy of cobalt and chromium was first 
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made in a high-frequency furnace. Using this as 


stock material, a series of runs was then made with 
various iron additions to determine by chemical 
analyses the losses met with in melting the same 
weight of material each time In order to compen- 
sate for these losses, desired amounts of the metal 
powders were added to the remelts before casting 
Before pouring, the melts were deoxidized with 
definite amounts of caleium-silicon alloy. The cast 
rods In. diam were then heat-treated 
Table 1 gives information about the heat 
ments and the chemical compositions of the cobalt- 
The chemical 


treat- 


iron-chromium samples investigated 
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ends of the samples. The iron contents varied from 
36.2 ? percent, and the chromium contents 
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from 8.5, to 9.8; percent by weight 


3. Apparatus 


The micrometric  thermal-expansion 
described by Hidnert and Souder [6] was used for 
the determinations of linear thermal expansion of 
the Rach 
was about 300 mm long 


apparatus 
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4. Results and Discussion 


Determinations of linear thermal expansion of the 
samples of cobalt-iron-chromium alloys were made 
at various temperatures between —65° and + 950° C 
and of the effeets of various heat treatments (from 

196 1,000° ©) on their thermal expansion. 
The results are summarized tn 3 tables and 8 figures. 


to 


Table 1 gives the coefficients of expansion that 
were determined for various temperature ranges 
from —60 300°C. The phase of the samples 
for each test and the residual change in length after 
each test are also given in table 1. 


Lo 


Figure 1 shows expansion curves of nine samples 
of cobalt-iron-chromium alloys between and 

or © These samples are arranged, from = the 
bottom to the Lop, in increasing values of their coef- 
ficients of expansion for the range from 20° to 
-60° CC. For all of the samples that were investi- 
gated between 65° and 60° C the observations 
were made in the following order: (a) Heating and 
cooling between room temperature and 60° © 
(test 1); (b) cooling and heating between room tem- 
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perature and 65° © (test 2); (¢) again 
and cooling between room termperature and +-60° C 
(test Several of the samples in figure 1 have 
coefficients of expansion that are very nearly zero 
for the range from 65° to +20° C., 

higure hows the expansion of a low-expansion 
alloy over the range from 65° io +300° C. While 
the average coefficient of expansion of this sample 
is less than 1 10-°/deg C over the range from —65° 
to + 60° C, it rapidly increases to 15.0 107*/deg C 
over the range from 200° to 300° C. It is seen that 
the low coefficients of expansion for such alloys 
generally apply for temperatures between 65° 
and +-60° C, 

Figure 3 indicates graphically the relations be- 
tween chemical composition and the average coef- 
ficients of expansion of the samples in the first tests 
for the range from 20° to 60° C. All of the coef- 
ficients of expansion apply for the samples in the 
original annealed conditions (y phase). The curve 
in this figure represents an isodil * equal to 1.0X 


heating 


> 
oO). 


a 


3’ The word “isodil’’ was derived in 1931 by Hidnert, from ‘iso’ (Greek isos, 
meaning equal) and from the first three letters of ‘dilatation’, and denotes a 
curve representing compositions having equal expansion, 
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TEMPERATURE, °C 


Linear thermal expansion of nine annealed cohalt- 
iron-chromium alloys (Fe 36.2, to 36.92, Cr 9.09 to 9.8:% 
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The initial observation for each alloy was taken at about 20° C and is plotted 


on the zero ordinate 
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Linear thermal expansion of annealed coball-tron- 
chromium alloy (Fe 36.58, Cr 9.06% 


. Heating; @, cooling;—, initial observation | 
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FIGURE 3. 
composition (percentage by weight) on the coefficients of linear 


er pansion (in millionths per degree by of annealed cobhalt- 


iron-chromium alloys for the range 20° to 60° C. 


10-° deg C, and appears to approximate part of an 
ellipse as shown by Masumoto [1]. The results 
indicate that the low-expanding alloys lie within 
this elliptical-shaped isodil and that the alloys out- 
side of it have coefficients that generally increase 
with increasing distance from it. These results 
indicate that close chemical control is necessary in 
the preparation of these low-expanding alloys. 

The samples with chromium contents between 8.5, 
to 9.1; percent have very low coefficients of expansion 
at room temperature but y--a@ transformations 
started between —10° and —S80° C. Ary tempera- 
tures for five samples may be found in table 1 (see 
also footnote q). These values indicate that for the 
samples investigated the Ar; temperature tends to 
decrease as the chromium content increases. No 
ya transformation was found on cooling a sample 
(see table 1, footnote r) containing 9.5; percent 
chromium, to as low as 196° ©. Masumoto [1] 
reported that his low-expanding allovs showed no 
ya transformations even when cooled to liquid 
hvdrogen temperature | 253° C), but four of his 
allovs expanded slightly (less than 0.02 percent) on 
cooling to low temperatures. 

Figure 4 shows the y--a@ transformation § that 
occurred on cooling sample 1762 to —67° C. This 
transformation was not completed after reaching 
this low temperature. The expansion on cooling 
caused by the y->e transformation proceeded with 
discrete rapid increases in length. These discrete 
increases in length were accompanied by audible 
clicks and slight increases in temperature of the 
sample. Upon heating the sample from —67° C. it 
was found that the rate of expansion was larger than 
it had been for the sample in its original annealed 
condition. It was found that for temperatures 
below 200° C the rate of expansion increases with 
increasing amounts of @ phase. 
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Ficure 4. 


Linear thermal erpansion curve of annealed cobalt- 
1ron-¢ hromiu m y 


alloy ke a. Cr 8.5,° 
transformation on cooling 


, Heating; @, cooling 


The ye transformation in sample 1766A (cut 
from same rod as sample 1766) proceeded 
thermally when its temperature was maintained at 
—70° C for 3 days. On rapid cooling from room 
temperature to —70° C, the sample expanded 0.12 
percent. Its length at —70° C increased with time 
at a diminishing rate. At the end of 3 days, the 
additional increase at —70° C was 0.06 percent. 

Figure 5 also shows an incomplete Ya transfor- 
mation that occurred in sample 1758 on cooling to 

64°C. Although the sample was cooled to —S80° 
C before the fifth test, the rate of expansion during 
the fifth test was practically the same as the rate 
during the third and fourth tests. This indicates 
that littl or no additional ya transformation 
occurred during cooling between —64° and —s0° C. 
The rate of expansion during the sixth test after 
cooling to 196° C. indieates that an additional 
ya transformation occurred in the sample. After 
the sixth test the sample was given a_ prolonged 
annealing at high temperatures (see table 1, footnote 
k) and then cooled slowly to room temperature. 
The smaller rate of expansion in the seventh test : 
compared to the rate in the sixth test of the sample, 
indicates that the annealing treatment resulted in a 
partial transformation of the atothe ycondition. Ref- 
erence should also be made to the coefficients of ex- 
pansion (between 20 ‘) in table 1 for a 


30.9 0) showing ya 


tial observation 


ISO- 


is 


and 60° C) 
better comprehension of the comparisons between 
the tests indicated in figure 5. 

Figure 6 shows the results obtained on six tests of 
sample 1760. The first three tests were made with 
the sample in the original y condition. The third 
test did not show any irregularities to 700°C. After 
the third test, the sample was cooled rapidly to 
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igure 5. Linear thermal expansion of annealed cobalt-iron- 


chromium alloy (Ke 36.9, Cr 8.86%). 


Combined tests 1 to 3 show y—a@ transformation on cooling; before fifth test, 
sample was cooled to —S0° C; before sixth test, sample was cooled to —196° C; 
before seventh test, sample was hydrogen annealed at 850° C for 10 hr and at 
1,000° C for 1 hr, and furnace cooled in 22 hr 

, Heating; @, cooling; —, initial observation 
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Figure 6 Linear thermal expansion of annealed cobalt-iron- 


chromium alloy (Ke 37.0;, Cr 8.879%). 
Before fourth test, sample was cooled to —196° C, 


, Heating Se. cooling 


190° C and held at this temperature several hours. 
A comparison of the expansion curves on heating 


in the third and fourth tests indicates that the 
cooling to — 190° C transformed a part or all of the 
to a An examination of the data in this paper 


indicates that the ye transformation in this sample 
was the most complete of any that were investigated 
in this work. On heating the sample in the fourth 
test, an a—»y transformation‘ started at about 600° C. 


‘Confirmed with other samples, 





Coefficients of linear thermal expansion of cobalt-iron-chromium alloy 


Fe 37.04, Cr 8.8; %) 


sample 1760; 


Average coefficients of expansion per degree C 


TABLE 2 
lest * Phase 
20° to 400° CoN 
3H Y 10. 4X10 11. { 
iC y 
4Hit a 11.7 12 
iC y+a 12.6 13 
5H >a 13 
H indicates that the coefficients of expansion were obtained on heating and 
Before this test, cooled to —190° C 


The ay transformation was completed during 
heating to 950° C in the fifth test, as indicated by 
the low expansion in the sixth test (compare with 
results of tests 1 to 3 for the y condition). 

Table 2 gives the average coefficients of expansion 
of sample 1760 for various temperature ranges from 
20° to 800° C. Figure 7 shows the coefficients of 
expansion versus temperature for the y and @ phases 
of this sample. The marked change in slope of the 
curve for the sample in the y condition at about 
240° C is typical for these alloys and is apparently 
associated with the magnetic transformation at 
the Curie point... For various mixtures of @ and y 
in this alloy, the corresponding curves for the coeffi- 
cients of expansion versus temperature generally lie 
between the curves shown for a@ and y (e. g., table 1, 
sample 1760). 

Figure 8 reproduced from Koster [7] indicates 
that the boundary between a@ and y phases is a steeply 
ascending curve up to 600° C in the region of 9 
percent of chromium. The conclusion, deduced 
from the data reported in table 1 of the present 
investigation, that the Ar; temperatures of the alloys 
decrease as the chromium content increases, confirms 
this equilibrium diagram. The Ac, temperatures 
obtained on three samples at about 600° C are also 
in agreement with this diagram. An examination 
of this figure in conjunction with data in table 1 
also indicates that @ phase should not be found in 
these alloys containing more than 9.6 percent 
chromium. 

‘The authors found that a specimen cut from sample 1757 in the y condition 
lost its magnetization at about 240° C Masumoto [1] indicated that the mag 


netic transformation point of low-expanding cobalt-iron-chromium alloys is tn 
the vicinity of 250° C 
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Ficure 7. Variation of coefficients of expansion versus tem- 
perature for y and a phases of a cobalt-iron-chromium alloy 
(Fe 37.0;, Cr 8.87%). 
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50° C 20° to 600° C 20° to 700° C20" to 800° C 
<10 13. 0«10- 13.9 10~ 
14.0 
15. 3X10 
14,2 14.9 15.3 


) indicates that they were obtained on cooling 


The coefficients of expansion in tables 1 and 2 
indicate the presence of a and y phases in the samples 
due to various heat treatments. Another indication 
of their presence is the residual change in length 
at 20° C. Table 3 gives the changes in length at 
20° C from the y condition of five samples after 
various heat treatments. The amount of @ phase 
present in a sample appears to be approximately 
proportional to its increase in length. 


3 Changes in length caused by transformations in 5 


TABLE 3. 
sam ple xs of cobalt-iron- hromium allo {8 


Cumulative changes * in length (at 20° C) after 
successive heat treatments 


Heat treatment 


1757 1758 170 THHA t L7SSA 

As received (+ condi- 

tion 0. 000 0. O00 0. OOO 0. 000 0. OOO 
Cooled to 65° C +. 170 +. 146 
Cooled to soe Cc 173 150 
Cooled to woe C +. 965 
Cooled to —196° C 179 23 +. 243 +. 182 
Heated to +700° C O67 070 
Cooled to 196° C O7s O70 
Heated to 1,000° C O69 OOS 006 OOl 
Held at —70° C for 

$} days 230 
Cooled to —196° C . 232 tT. 182 


*» Where values are not indicated, the heat treatments were not given, 
b Cut from same rod as sample 1766 
¢ Cut from same rod as sample 1783. 
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Two sections of ternary equilibrium diagram of 
alloys for 40 and 30 percent iron 


FIGURE 8. 
cobalt-iron-chromium 
(Késter [7]). 
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Microstructure (75 
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FIGURE 9, 


Etched with hydrochloric and picrie acids in aleohol [10] 


A metallographic examination ® was made of some 
of the samples in the y and a+y conditions. Typical 
structures are shown in the photomicrographs of 
figure 9. This figure shows that the annealed 
sample (fig. 9A, y) has large grains and dendrites 
whereas the sample that had been cooled to —196° C 
(fig. 9B, a+y) has a needle-like structure. The 
needle-like structure formed during the transforma- 
tion on cooling and the kinetics of the transformation 
indicated by the dilatometric experiments, shows 
that this transformation is of the martensitic type 
[8,9]. 


5. Summary 


1. The linear thermal expansion and phase trans- 
formations of some cobalt-iron-chromium alloys 
were investigated at various termperatures between 
—65° and + 950° C. 

2. Annealed alloys having chemical compositions 
of about 54 percent cobalt, 36.6 percent iron, and 
8.9 percent chromium have coefficients of expansion 
less than 1X10~-"/deg C for the range from 20° to 
60° C. Close control of chemical composition is a 
prerequisite for success in the production of low- 
expanding cobalt-iron-chromium alloys. 

3. For temperatures above 60° C, the coefficients 


”. 
of expansion of the alloys increase to about 16x 107°/ 
deg C for the range from 250° to 300° C 

4. Some low-expanding alloys were found to have 
a yea transformation on cooling to low tempera- 
tures. This transformation of the martensitic 
type in that it produced a needle-like structure and 


A 


Is 


6 By D. B. Ballard of the Metallurgy Division of the National Bureau of 
Standards 


of two cobalt-iron-chromium alloys. 


Sample 1766, after annealing at 1,000 


| 
| 
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C; B, sample 1757, after cooling to 


proceeded with discrete rapid increases in length, ac- 
companied by audible clicks and slight increases in 
temperature. 

5. The coefficients of expansion below 200° C in- 
crease with increasing amounts of a phase present 
in the alloys. 

6. Samples containing @ started to transform to 
y at about 600° C, 
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Vapor Pressures of the Methanes 
George T. Armstrong, F. G. Brickwedde, and R. B. Scott 


Experimental measurements, not previously reported in detail, of the vapor pressure of 
CH, and the differences in vapor pressures of CH, and the deuteromethanes below the normal 


boiling point are presented. 


of methane and the deuteromethanes is give 


n. 


A critical review of the literature relating to the vapor pressures 


On the basis of the experimental measure- 


ments and the previous literature the following best values are assigned to fixed points on the 


vapor pressure curves. The triple points are: CH,, 87.50 mm, 90.66°; CH;D, 84.5 mm, 
90.41°; CH,D,, 82.0 mm, 90.17°; CHDs, 80.2 mm. 89.96°; CD,, 79.1 mm, 89.79° K. The 
boiling point of CH, is 111.67° K; the critical point 45.6 atmospheres and 190.6° KK. Equa- 


tions and a table are given for the 


vapor pressure of CH, in the solid region, the liquid region 


below the normal boiling point and the liquid region above the normal boiling point. 
yr ratios of the vapor pressures of the methanes are expressed by equations of the 


form T logy Pp/Puy=A—B/T. The constants A and B have the following values, respec- 
tively: in the solid range CH;D, 1.260, 110.2; CH.D., 2.694, 222.2; CHDs, 4. a 351.7; 
CD,, 5.529, 410.5; in the liquid range, CH;D, 1.328, 129.5; CH,.D»:, 2.671, 245.4; CHDs, 
3.969, 343.8; CD,, 5.159, 421.1. The vapor-pressure ratios of the isotopic me a forma 
nearly geometric progression with increasing deuterium substitution. Deviations of their 


mixtures from ideal solutions are very small. 


1. A Review of the Vapor Pressure of Meth- 


ane and the Deuteromethanes 
1.1. Methane 


vapor wire - solid and liquid methane 
has been determined in various temperature ranges 
by many sac aera png n addition, there have 
been several summaries of the published data. S. F. 
Pickering [1]? reviewed the literature relating to the 
critical constants; Copson and Frolich [2] summarized 
the vapor-pressure data available in 1929; and Stull 
summarized the data in 1950 [3]. A series of best 
oe appeared in the International Critical Tables 
[4]; Egloff [5] and Timmermans |6] surveyed the 
literature in preparation for their reference books of 
physical constants of hydrocarbons; and a table of 
selected values for the vapor pressure below the 
normal boiling point was given by Rossini and co- 
workers in collaboration with the American Petro- 
leum Institute [7]. Other, less critical, surveys have 
also been made. An important deficiency remaining 
after the above summaries that the excellent 


The v 


Is 


API table of vapor pressures in the low-pressure — 


range is based to a considerable extent on data not 
generally available. 

In an attempt to remedy this deficiency this 
view has been prepared in which the data from 
various sources have been compared including those 
described in section 2, previously unpublished, and 
available to none of the previous reviewers except 
perhaps Rossini. An attempt has been made to 
find equations which fit well the best available data 


and at the same time pass through selected best 
values for the triple point, normal boiling point, 
and critical point. 
This research was supported in part"by the’Atomic Energy Commission 
Figures in brackets indicate the literature references at the end of this paper 
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The vapor pressure of solid methane has been 
measured by Hunter [8], Henning and Stock [9], 
Karwat [10], Freeth and Verschoyle [11], Tichner 
and Lossing [12] and by Brickwedde and Scott of 
the National Bureau of Standards (section 2 of this 
paper). Determinations of the triple point as listed 
in table 1 provide a fixed point through which an 
equation for the data should pass. Of the reported 


determinations, the extensive data from Clusius’ 
laboratory appear to be best represented by the 
values 90.67 +0.03° K and 87.4 +0.1 mm, which 


were determined from a study of the melting tem- 


perature as a function of pressure. The values 
90.660 +-0.005° K and 87.60 +0.10 mm. were found 


by Brickwedde and Scott with a carefully purified 
sample. The best values of the triple point tempera- 
ture and pressure appear to be 90.66 +0.01° K and 
87.50 +0.10 mm. 


Taste 1. Triple point of methane 
Observer Date Temperature Pressure 
Kk mm He 
Olszewski [13] ISS 87.4 sO) | 
Hunter [Ss] 1u0t WO. 5 to V1.2 y2 | 
Crommelin [14] 1912 OO). O1 70 
Henning and Stock [9] 1921 YU. 0 78.3 | 
Eucken and Karwat [15] 1924 WO. 6 x7 | 
Clusius [16] 1929 Ww. ¢ | 
Wiebe and Brevoort [17] 1930 Ww 
Freeth and Verschoyle [11] 1931 W). 70 b 87. 5240.1 
Fischer and Klemm [18] 1935 Ww. 7 ASN 
Kruis, Popp, and Clusius [19] 1937 00.6 87. 5 
Clusius, Popp, and Frank [20] 1937 WO. 62 87.5 
Frank and Clusius [21] 1937 00. 6 87. 440.05 
Brickwedde and Scott (this re 
port 1037 00. 660-40. 005 87. 604-0. 10 
Clusius and Wiegand [22] 1940 00. 674-0. 03 87.440. 1 
Stavely and Gupta [23] 194y 87. 7540.07 | 
Selected values 00. 66 R7. 50 


» Corrected by 0.1° K for change of temperature scale. 
Average of two values 








Of the various equations that have been proposed 
the Antoine type equation given by Rossini et al., 
appears to give the most satisfactory fit to the data 
from Various sources 

This equation presents a compromise between the 
NBS data and those of Freeth and Verschoy le and of 
Karwat in the region from SO° kK to the triple point, 
at which point it represents a good average of the 
pressures found by Brickwedde and Seott and b 
Clusius and Wiegand The equation is in excellent 
agreement with the recent smoothed data of Tichner 


by 


and Lossing down to 0.01 mm Below this it devi- 
ates by as much as 10 percent from thei data 
When eq (la) given by Rossini 
log)o/? (mm) =7.69540— 532.20/( 7+ 1.84 La 


») 


is solved simultaneously with eq for the liquid 
state vapor pressure, these two equations give values 
90.67° IKK and 87.7 mm for the triple point tem- 


perature and pressure. In the original equations by 


of 


Rossini temperature Is viven In degrees +. and for 
the conversion to absolute temperature (7°) the ice 
point is taken as 273.16° K. As the value of the 


pressure is higher than the reported values in general, 
a very small adjustment was made in the equation 
for the solid state vapor pressure to bring the triple 
point pressure 87.50 mm. The corresponding 
temperature ts 90.66° in agreement with the value 
selected above The revised equation for the vapor 


to 


The adjustment is so slight that it changes the value 
of P at any given temperature by a maximum of 
0.02 mm, which is well within the limits of error of 
most of the data. Besides the achievement of con- 
sistency at the triple point of CH, this adjustment 
appears also to be justified by a greater consistency 
in the vapor of the deuteromethanes at 
their triple points, as is shown later. The individual 
vapor pressure measurements of various investigators 
have been compared with eq Ib) and the deviations 
are plotted in figure 1. 


pressures 


The vapor pressure of liquid methane in the range 
below its normal boiling point has been represented 
in the published literature 
[4, 24, °5]) most frequently by tables based on the 
relatively few of Henning 
and Stock [9]. Tavlor, and Smith 
[26] lving 1 or 2 percent higher, but of equally good 
internal consistency 


for example: references 


reported measurements 
Data of Keves, 
in this region, have generally 
Besides these two sources of data 
of Hunter [8] 
A few 
were obtained by Cragoe [27] in this region but were 
never published. To these may be added the data 
in 2. which internal 
precision and cover the range more thoroughly than 


been disregarded 


there are some carly measurements 


which appear to be of low accuracy values 


described section are of high 


previous studies. Equation (2) published by Rossini 


- é et al., appears to be based largely on the latter data. 
pressure of solid methane is : 
log, .J?(mm) =7.69540 — 532.20/(7+- 1.842). (1b log,,/?(mm) = 6.61184—389.93 /( 7 —7.16 2) 
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The normal boiling point calculated from this 
equation is 111.67° K.  Beeause this value is only 
0.001° from the value found by Brickwedde and 
Scott and within 0.01° of the mean of the values 
reported by Henning and Stock, by Cragoe, and by 
Keves, Taylor, and Smith, it appears to be the best 
value, with an estimated uncertainty of +0.01° K. 
The uncertainties listed here and later are estimates 
of the limits within which the true value has a 50 
percent probability of lying. Observed values are 
shown for comparison in table 2. The deviations of 
the measured data from eq (2) are shown in figure 1. 
It will be seen that the equation agrees closely with 
the data of Brickwedde and Seott and les near that 
of Henning and Stock at the center of the range. It 
may be noted here that in addition to eq (1) and (2 
other equations of different functional form can be 
made to fit the data with nearly as vreat precision, 
as for instance those given in section 2. 


PARBLI 2 Vo mal ho ng pon methane 
O ) Ie ‘ 
A 
‘> ew sa ] 
Hurt s nT iW 
He ist nr 111.79 
Keyes, 1 > y22 111 
~. 5. 4 e [2 ng 111. 5S 
Brick we t - ) 111. 669-0. 010 
Selected I11.t 


In the pressure region above one atmosphere there 
are extensive but widely divergent data. Certain very 
early data such as those of Wroblewski [28] and of 
Olszewski [13] may be disregarded as lacking suf- 
ficient: precision when compared with later data. 
There remain the data of Cardoso [29], Cragoe [27], 
Keves, Taylor, and Smith [26], Eucken and Berger 
30], Volova [31], and Bloomer and Parent [32]. 
These also differ widely among themselves. The 
paper of Eucken and Berger, for example, does not 
include any description of apparatus or technique 
and is reported in such a way as to leave the impres- 
sion that the data are of only moderate accuracy. 
The seatter of the points is quite bad particularly 
at the lower temperatures. The data of Cardoso 
covering a limited temperature range are of a precision 
comparable with that of Keves and of Cragoe, but 
differ from theirs by about 2 percent. Keves and 
that an error of 0.56° in 
Cardoso’s temperature scale would account for the 
discrepancy, but as was pointed out by Pickering 
the line of argument leads to a temperature change, 
which if properly applied, widens the error instead 
of closing it. 

The correction sugeested by Keves involved the 
use by Cardoso of the melting point of toluene as a 
calibration point for his thermometers, and a change 
in the aecepted value of the melting point which 
occurred after Cardoso’s measurements were pub- 
lished. Cardoso used —94.5° C for this temperature. 
The thermometers used by Keves were also used to 
measure the melting point of toluene [33], —95.7° C 
being found. It might be pointed out that if this 


coworkers suggested 


4l 


argument were extended to the present and the best 
value accepted for the melting point of toluene 

94.991° C as found by Rossini et al. [21], then it 
appears (if the purity of the toluene and the tech 
nique of determining melting points are assumed to 


be sufficiently good in the early investigations) that 
the true temperature scale must’ be somewhere 
between that used by Cardoso and that used by 
Keves. The assumption concerning the purity of 


the toluene is of doubtful validity. Assuming i 
to be true for the purpose of discussion, Cardoso’s 
pressure measurements belong to lower true temper 
atures than he recorded. This results in raising the 
vapor pressure values for a particular true tempera- 
ture. Similarly, Keves’ pressure measurements be- 
long to higher true temperatures than he recorded, 
or the vapor pressure for a particular true tem- 
perature should be lower. This would widen the 
gap existing between the data by an amount cor- 
responding to 1.2° at the eritical point. Adding this 
to the existing difference of about 0.5°, the difference 
in the data come to about 1.7° C near the critical 
point, or for a given temperature to about 5 percent 
in pressure. Examination of the several melting 
points recorded by Keves and coworkers [33] sug- 
gests the possibility of a consistent error of 0.6° to 
O.8° in his temperature seale in the region from 

112° to —95° C. 

The observed critical point pressure should be 
independent of temperature scales. Here we find 
that Cardoso measured a value 0.46 atm. lower than 
that of Keves. The difference ean be due either to 
a difference in the purity of the CH, used, or to eali- 
bration errors in the pressure scale. In this respect 
the calibrated steel tape used by Keves for measure- 
ment of mercury column height seems definitely 
superior to the volume technique used by 
Cardoso. 

In addition to the published data for methane, 
some vapor pressure measurements were made by 
Cragoe [27] in 1919, which were never published. 
These show an undesirable amount of scatter, but 
in general tend to fit the data of Keyes, Taylor, and 
Smith. Cragoe also tabulated values above 1 atm. 
for the International Critical Tables which are based 
primarily on his own data and that of Keyes. 

Volova [31] in a study of phase equilibria in 
methane-ethylene mixtures reported values for the 
vapor pressure of CH, for which he claimed an 
accuracy of 0.2° in temperature. The values for 
pressure which he gives lie within this limit of eq (3) 
and much closer to the data of Keyes and of Cragoe 
than to any other data. 


gas 


TABLE 3. Critical point of methane 


Observer Date Temperature Pressure 
kK atm 
Dewar [34 Iss4 173.7 17.6 
Wroblewski [2s Iss4 190. 7 MLS 
Olszewski [13 ISS5 1.4 4.0 
Cardoso [29] V5 100. 31 5.00 
Keyes, Taylor, and Smith [26] 1922 191,06 3. 06 
Bennewitz and Andreev [35] 1920 100. 7 45.7 
Bloomer and Parent [32] 1952 190, 55 45.47 
Selected values 100. 6 5.6 








TaBLe 4. Vapor pressure of methane 
T°’ K logy P P(mm)* T°K | logwP | P(mm)s 
Solid Liquid 
aD | >7. 62387 0. 0042 121 3. 18910 1545. 6 
52 7. 81093 0065 122 3.21914 1656. 3 
53 7. 99116 0098 123 3. 24865 1772.8 
“4 8. 1404 0146 124 3. 27766 1895. 2 
55 8. 33261 0215 125 3. 30618 2023. 9 
56 8. 49448 0313 126 3. 33422 2158.9 
57 8. 65084 O445 127 3. 36179 300. 3 
58 &. SO199 06384 128 2. SSS892 1448. 6 
59 8. 94816 OSR7 129 3. 41561 2603. 8 
60 9. 08960 1229 130 3. 44187 TH6. 1 
6l 9. 22655 1685 131 3. 46771 2935 
62 9. 35920 2287 132 3. 49316 3112.9 
63 9. 48776 3074 133 3. 51821 3297.6 
4 9.61242 410 134 3. 54287 3490. 3 
65 9. 73335 ‘41 135 3. 56717 S691. 3 
66 9. 85071 709 136 3. 59111 3900. 4 
67 9. 96466 922 137 3. 61469 4118.0 
68 0. O7534 1. 189 138 3. 63793 4344.4 
} 69 18200 1. 524 139 3. 66083 4579.6 
70 28747 1. 939 140 3. 68341 $824.0 
71 38917 2. 450 141 3. 70567 5077.8 
72 $8511 3. O77 142 3. 72762 5341.0 
73 SS441 3. S41 143 3. 74927 5614.0 
74 67817 4. 706 144 3. 77063 5897. 0 
75 76950 5. 88 145 3.79170 6190. 1 
76 S547 7.22 146 3. 81248 (493. 5 
77 M4519 S81 147 3. 83300 6807.7 
78 1. 02973 10. 71 14s 3. 85325 7132. 7 
7¥ 1.11218 12. 95 149 3. 87325 7469.8 
SO) 1. 19262 15. S58 1M) 3. SY2UR 7816.0 
Sl 1. 27112 18. 67 151 3. 91247 8174.7 
S2 1. 34774 22. 27 152 3. 93173 8545.3 
AS 1. 42256 26. 46 153 3.95074 S028, 8 
4 1. 49563 31. 31 14 3. WHY5S Y322. 4 
SS 1. 56702 36. 90 155 3. GRSDO Q729. 5 
SH 1. 63679 $3.33 156 4. O42 10149 
s7 1. 70499 “70 157 +. 02456 10582 
ss 1. 77166 59. 11 158 4.04248 11028 
Su 1. 83687 68. 69 159 4 OH019 11487 
uO 1. 90086 79. 5S 160 1.07771 11959 
161 4+. 09502 12446 
Liquid 162 4.11215 12046 
163 4. 12909 13461 
It 4. 14585 13091 
ul 1. went 41. 40 165 4. 16242 14535 
v2 2. O1578 103. 70 166 4. 17882 15004 
03 2. 06932 117. 31 167 1.19506 L670 
bt | 2. 12163 132. 32 168 4.21112 16200) 
U5 2.17275 148.85 16a 4. 22702 1HS66 
o6 2. 22271 167. 00 170 4. 24276 17489 
v7 2. 27157 1S#, SS 
Us 2.31935 2s. 62 171 4. 25834 18128 
ay 2. 36609 232. 32 172 4. 27376 INTS2 
100 2. 41182 258. 12 173 4. 28993 19455 
174 1.30415 20144 
101 2. 45658 Zs. 14 175 1.31913 20851 
102 2. S039 316. 51 176 4. 333906 PTA 
103 2. 54329 349. 3S 177 +. 34806 22318 
104 2. 58530 384.85 178 4. 36321 23078 
10 2. $2046 $23. 12 179 4. 37704 23859 
106 2. HH678 4. 28 Iso) 4. 39192 2466 
107 2. 70629 508. 50 
108 2. 74502 555. 93 18] 4. 40608 25473 
109 2. 78200 606. 73 182 4. 42011 26310 
110 2. 82022 661. 01 183 4. 43402 27161 
| 1N4 4. 44780 28041 
lil 2. 85674 719. 02 ISS 4. 46145 28037 
112 2. 89261 780. 92 186 +. 47499 29853 
} 113 2. 92814 S47. 50 187 4. 48842 s9791 
114 2. 96204 918. 20 ISS 4. 50172 31748 
115 2. 99707 4935. JS 1SY 4.51491 32727 
116 3. 0305 1072.9 140 4. 52799 43728 
117 3. 06340 1157.3 
118 3. 09572 1246. 6 
119 3. 12741 1340.9 
120 3. 15853 1446.6 


* The final figures in the values of P(mm) in general are not significant 
> Logarithms have been increased by 10 where necessary to avoid negative 
characteristics 


The very recent work of Bloomer and Parent [32] 
provides good confirmation of the data of Keyes in 
the region between the boiling point and about 180° 
K. Above that it tends to confirm the suggestion 


| 


those of Keyes and Cardoso. Because it was re- 
ceived after the selection of the values listed in table 
4, it did not receive any weight in the establishment 
of these values, though it appears to offer a distinct 
improvement in the high pressure region. As was 
noted by Timmermans the critical constants of 
methane have been extremely uncertain considering 
the importance of this substance. The values 190.6° 
K and 45.2 atm satisfy eq (3). However, a better 
value for the pressure on the basis of recent measure- 
ments would be 45.6 atm. Table 3 lists measured 
values of the critical pressure and temperature for 
comparison. For purposes of interpolation the 
vapor pressure equation of Keyes’ is unsatisfactory 
in that it does not join well with the recent low 
pressure data. Correction terms were added to the 
Keyes equation which bring the equation close to the 
values for pressure given in the International Critical 
Tables, and at the same time cause the equation to 
pass through the normal boiling point, 111.67°. 
Equation (3), derived in this way, was used as a 
basis for calculating deviations of experimental data 
plotted in figure 1 in the range above 111.67°. 


logy/?(mm) = 10.6863 1 —595.546/ T—0.0348066 T 
+-0.00013338 T? — 1.7869 X 107° T! (3) 
Errors larger than 1 percent are unlikely in view of 
the values reported by Volova, by Cragoe, and by 
Bloomer and Parent. 

As a summary of the vapor pressure of CH,, 
values have been calculated at selected temperatures 
with the use of eq (1b), (2), or (3) as appropriate, 
and these values are shown in table 4. 


1.2. Deuteromethanes 


Aside from the studies made by Brickwedde and 
Scott at the Bureau, which have been briefly reported 
in a preliminary form [36] and which are deseribed in 
detail in section 2, the vapor pressure data on the 
deuteromethanes are rather sparse. A boiling point 
determination by Stedman [37] on a mixture of CH, 
and CH,D indicated that the vapor pressure of 
CH,D is higher and that it might boil 0.5° lower than 
CH,. A number of measurements were made of the 
triple point temperature and pressure of CH,D and 
CD, in Clusius’ laboratory [19, 20, 22, 38] and in his 
laboratory a series of determinations of the vapor 
pressure of CH,D was made over a very narrow 
range [21]. In table 5 is shown a comparison of the 
reported triple point temperatures and pressures. 
There is very good agreement between the values 
found at the Bureau, as calculated on the basis of 
a recent analysis of the samples, and those reported 
by Clusius and Wiegand which appear to be the most 
reliable results from Clusius’ laboratory. In figure 
2 are shown the vapor pressure data on the deutero- 
methanes, plotted as the function 7 log jl?p/Px 
where Pp is the vapor pressure of the deuterium 


that the true temperatures lie somewhere between | compound and Py the vapor pressure of CH, at the 
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TaBLe 5. Triple point temperatures and pressures of the 
deuteromethanes 
CH;D CHD» CHD CDs 
Observer 
AT*|AP*|AT|AP|AT|AP|AT| AP 
| = 
Kruis, Popp, and Clvsirs 
{19 1. 46 8.6 
Clusius, Popp, and Frank 
[20] 0.C8& 3.1 1. 36 8.6 
Clusis and Wiegan? [22] 24 | 3.1 SS 8.6 
Clusius and Pop [38] 24. 3.0 SS &. 6 
Brickwelde and Scott 
(section 2 254 2.08 0.492 5.470.701 7.31 BOS 38 
Selected best values 25 | 3.0 499 5.5 70 |7.3 87 4 
*The trimle point temperature of the deuteromethane is 90.66—AT7 °K, and 
the triple point pressure of the deuteromethane is 87.50—AP mm 
TABLE 6. Constants for the equations T logyPp/Pa==A— B/T 
Solid range Liquid range 
Pure substance e« (mm 
i B A B 
CH.D 1. 260 110.2 1. 328 120. 5 +0. 015 
CHD» 2. 604 222. 2 2. 671 245.4 O17 
CHD 4.452 s51.7 $. 409 343.8 OLS 
CDs 5. 429 410 159 $21.1 aoy 


same temperature. The experimental NBS values 
are shown as circles; the area covered by the data of 
Clusius and Wiegand on CH,D is shown by a dotted 
rectangle. The spread of their data overlaps the 
NBS data, though the average is higher. The large 
spread represented by their data is probably a result 
of measurement of absolute vapor pressures only, 
with an accuracy of 0.1 mm. A considerably more 
accurate measurement of the differences between the 
methanes was achieved at the NBS by the use of 
differential oil manometers. 

Best values for the triple point temperatures and 
pressures were selected by comparison of the values 
given in table 5. For values of the vapor pressure 
at other temperatures, in the absence of other data, 
constants for eq (4) have been derived as described 
in section 2 on the basis of the NBS data alone, 

T logwPp/Pa=A—! (4) 
using analysis 5 for the compositions of the samples. 
The constants for CD, are derived on the bais of 
measurements on sample 1 only. 

The constants “A and / derived in this manner are 
shown in table 6. These constants differ slightly 
from those previously reported [36] which were 
based on the use of the earlier less reliable analysis 4. 
The mean deviations of the experimental points 
from these equations are also shown, e being the 
deviation of AP» 4») In mm Hg. One source of 
uncertainty whose effect is not indicated by the 
deviations is the uncertainty of composition of the 
samples. The good agreement between the cal- 
culated values for the vapor pressure of CD,, based 
upon measurements made on two samples of widely 
differing composition, is an indication that the uncer- 
tainty due to lack of knowledge of the composition 
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FIGURE Vapor pressures of the pure deuteromethanes. 


of the CD, samples is very small. The greatest 
uncertainty in the vapor pressure difference 
AP cup cH, is due to the lower Jimit of detectability 
of CH, in the presence of CH,D by mass spectrom- 
eter, which is approximately two percent. The com- 
position induced uncertainties in the vapor pressure 
differences AP cx.p,-cay and AP cyp,-cu,y are some- 
what greater. While this uncertainty is difficult to 
assess, an estimate of an upper limit to it may be 
obtained by comparing the constants found as the 
results of several different attempts at analysis, and 
shown in table 10 of section 2. 


2. Experimental Determination of the Vapor 
Pressures of Methane and the Deutero- 
methanes 


The experimental work on the measurement of the 
vapor pressures was carried out in the years 1935 to 
1938 by F. G. Brickwedde and R. B. Scott. After 
the experiments were completed the samples of the 
deuteromethanes were retained in storage. Because 
new information for the interpretation of the ob- 
servations became available, a new examination of 








the samples was undertaken in 1952 and a reevalua- 


tion of the data was carried out by G. T. Armstrong 


2.1. Preparation of the Methane 


Six 4-liter sumples of methane were prepared as 
indicated by the following reactions 
9CO+2D,.=—=CD,+CO 
2) CCL+4C.H.D — 4Zn = CD, + 4C0,H-OZnCI 
CHCL + 3C.H.0D + 3Zn— CHD, + 3C,H.0ZnC! 
1) CH. 1L,+2C,.H.0D + 2Zn=CH.D.+ 2C.H-OZ4n!I 
5) CH.Mel~ D.O-— CH,.D~ DOMe!I 
Hb) CCT (CHO - 4Zn—CH 1(C.H-OZAnC! 


were prepared by Dr 
Sample 


All of these samples except | 
Kivoshi Morikawa at Princeton University 
s prepared by Weber at Columbia University 
The nonmethane impuri removed from these 
sumples by fractional distillation in a small low 
temperature still. For this work the still was sur 
rounded by a liquid oxygen bath boiling at about 81 
em He Two distillations made ‘h 
San] le Krom the original 4-liter sample a middle 
fraction of 1.5 to 2 liters was collected in the first dis 
tillation and 300 to 500 cm? in the second distillation, 
Before and after the collection of the middle 
“pure” fraction of the second distillation, 200 em 
fractions were collected The vapor 
these 200 cm? fractions at their boiling points did not 
differ from that of the middle fraction by as much as 
0.1 mm Hg and the difference was usually much 
smaller. This was regarded as evidence that im- 
purities whose vapor pressures were very different 
from that of methane were not present in significant 
amounts in the “pure” fraction. 

Another sample of CH, was obtained by the frac- 
tional distillation of some relatively pure methane 
purchased from the Linde Air Products Co. This 
sample, 1,900 ¢m*, was used in all the vapor pressure 
work. The vapor pressure of this sample differed by 
less than 0.1 mm Hg from the vapor pressures of the 
200 em® samples collected immediately before and 
after it during the distillation. The vapor pressure 
of this material differed from that of the synthetic 
CH, sample listed above by less than 0.1 mm He. 


1 wa 
ics were 


were on Cit 


or 


pressures of 


2.2. Apparatus 


The vapor pressure of CH, and the differences be- 
tween the vapor pressures of the | eavy methanes and 
CH, were measured by means of the vapor pressure 
apparatus shown in figure 3. The copper block F, 
containing three into which the methane 
samples were condensed, bangs in an evacuated 
vessel, D, which is surrounded by a bath of liquid air 


cavities 


in the Dewar flask, I. The insulating vacuum in 
vessel D is produced by a vacuum pump at A. The 
thin-walled copper nickel tubes, } connect” the 


cavities to the differential oil-manometer, J. One of 
the tubes also connects the absolute mereury 
manometer, K. The temperature of the block ts 
given by the small platinum resistance thermometer, 


to 
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FIGURE 3. Di fhe ential vapor pressure apparatu 
A, V, Vacuum pump conn B, filling tubs C, thin-wal 
ulbs D, evacuated ve l; E, ther i] eld; 1 pper block; G, platinun 
re tance therrmometer; Hl, constant temperature bith enclosure; I De ul 
flask: J, differential oil manometer; kK, absolute mercury manometer 
G. This thermometer is one of the standards which 


has been compared with the gas thermometer [39] be- 
tween 10° and 90° K. Above 90° K the Interna- 
tional Temperature Scale calibration was used. The 
temperature of the block is maintained uniform and 
constant by controlling the temperature of the 
thermal shield, E, which is equipped with electrical 
heaters and differential thermocouples. The tem- 
perature drift of the block was less than 0.0001 deg 
per minute. To guard against cold the 
tubes, C, tl ey were wound with an electrical heater 
from a point at the top of the shield to a point above 
the level of the bath. <A differential thermocouple, 
placed just below the lower end of the heater, was 
kept slightly warmer than the block. 

The lowest temperatures, down to 54 
tained by reducing the pressure in vessel H, contain- 
ing the liquid air bath, by means of a large vacuum 
pump connected at V. Hydrogen at a few milli- 
meters pressure was let into the insulating vacuum 
space while cooling the apparatus. It was found 
necessary, when working at the lower temperatures, 
to leave a little hydrogen in the insulating vacuum 
space to conduct away the heat flowing down the 
copper-nickel tubes. The various samples were 
introduced into the apparatus through the tubes B. 
Connections not shown in the diagram made pos- 
sible the thorough evacuation of all the apparatus, 
so that the samples were not contaminated, A mer- 
cury displacement pump was employed in admitting 
and removing the samples. 


spots on 


K. were ob- 


The vapor pressure readings were taken on mirror 
scales supported against the backs of the ma- 
nometers. These scales could be read to 0.1 mm, 
The absolute mercury manometer, K, gave the vapor 
pressure of CH,; the differential oil manometer J 
measured the vapor pressure difference, heavy 
methane minus CH,. The oil in J is a diffusion 
pump oil having a very low vapor pressure and is 
chemically stable when heated to drive out 
solved gases. As an extra precaution against con- 
tamination of the liquid methane, while in the vapor 
pressure apparatus, by oil vapor, traps surrounded 
by solid CO, were placed in the lines leading from 
the block to the manometers. 


dis- 


2.3. Vapor Pressure-Temperature Relations for CH, 


The vapor pressure of CH, was measured in the 
temperature range 54° to 112° K. If was found 
that the observations could be represented by the 
following equations 


to 90.66° K 
177.46/ 7-4 (5) 


CH, (solid): temperature range 54 
logio/?(mm Hg) 6.7838 
0.00469 7 


CH, (liquid): temperature range 90.66° to K 112 
logy/?(mm Hg) =7.55073 
183.22 7'— 0.0030686 7’. (6) 


The temperatures recorded on the International 
Temperature Scale were converted to degrees Kelvin 
using 7)=273.16° Kk. The deviations of the obser- 
vations from the equations are shown in figure 4. 
These deviations Than be considered as representing 








reliable in the region hear YO 
K since the thermometer, while in the 
apparatus, was calibrated at the point of 
OXVeen. Pure oxveen ‘was condensed into the y upor 
pressure apparatus and vapor 
readings were taken in the neighborhood of the OX\ 
gen boiling point. Later calibrations of the resist 
ance thermometer in the upparatus customarily used 
for oxygen boiling point calibrations on the [Inter 
national Temperature Seale, results) which 
agreed to 0.005 deg with the calibration in the 
methane apparatus. 

A reasonable value for 
peratures near 90° Kis +0.005°. At 111.6 
estimated uncertainty ts +0.010° KO which 
sponds to an error in vapor pressure of + 0.6 mm Hg. 
At temperatures considerably below 90° K the un- 
certainty increases to -+0.02°, but this is of no con 
sequence since dP/dT is so small that the over-all 
accuracy is limited by errors in the measurement of 
vapor The uncertainties in temperature 
given above do not take into account any error in 
the accepted value of the oxygen boiling pot, 
which is taken to be 90.19° K or 182.97° C Like- 
wise, above the oxvgen boiling pot no consideration 
was given the possible deviations of the International 
‘Temperature Seale from the thermodynamic seale. 


given should be very 
resistance 


boiling 


pressure resistance 


rave 


the uncertainty at tem 
K the 


COTTC- 


pressure. 


2.4. Vapor Pressure Difference (CH,,D,_ ,,—-CH,) 


The differences between the \ wpor pressures of the 
five samples of heavy methane and that of CH, were 
investigated in the range of temperatures 54° to 112° 
K. The data are represented in the graphs, figure 5. 


4 
the errors in the measurement of Vapor pressure since 
the precision of the resistance thermometer readings . - : 4 
° ‘ e T J 
was such that the expected resultant scattering would t 
: pnts A ee , 
be only 0.06 mm Hg at the boiling point of CH, and 4 ; 
much less at lower temperatures. Constant errors Ya ; 
. " ° ° / 4 / 
resulting from inaccuracy of the temperature scale, M4 ; 
: ry |} f+ + s 
should be considered separately. rhe temperatures - 
warencticon VA A 
Phese equations deviate no more than 0.1 mm from eq (1) and In Section f \ ‘ i : é 
1, although ditferent in for 4 
‘This oxvgen Was pre ired by the thermal decomposition of thoroughly out o/ . Y 
gassed potassium permanganate Its vanror pressure was compared with that ie | ont | 
ofasample prepared by the electrolytic dissociation of a carefully purified solution wit hinet— —— % — r 
of baritim hydroxide, Shepherd, Weaver, and Pickering, J. Research NBS 22 v 
01 (1939) RPLIS2. There was no nificant difference in the vapor pressures 4 ? 
f the two sampk kk than 0.1 mm Hg Z y e 
€ 
J > 
TRIPLE Jf , 6 
' 
POIN 4 
2 k : 
] , 4 A+ 8 e 
} P y 5 
¢ 
| 
> | yy a 
. P f a 
é - ha { pe H 
tal ’ 
| 
1 | os r “ + 
SX 
| | i 4 P ‘ 
i oa 
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At uw v "3 
| yw » a a 
| iy : i a pe 
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FIGURE 4, Deviation obs-calc) of CH, vapor pressure measure- PiGure 5, Vapor Pressures oj impure deuteromethane sample Ss. 
ments from eq )) and (6), | Ihe insert is an enlargement of the low pressure range; 
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The deviations of the observed points from the 
smooth curves average less than 0.01 mm Hg. It will 
be noted that there are large differences between the 
values obtained for the two samples of CD,. A 
logical explanation of this discrepancy is that none 
of the samples were pure, but consisted of mixtures 
of the several isotopic modifications of methane. 


2.5. Triple Point Pressures and Temperatures 


The triple point pressures of the methanes were 
measured in the following manner: The sample was 
condensed into a spherical pyrex bulb of about 10 
cm volume at the end of a tube 30 em long, 5 mm 
inside diameter. This tube communicated with a 
mercury manometer. The sample was frozen with 
liquid air, then the liquid air flask was lowered until 
the vapor pressure started rising slowly. When the 
triple point was reached, the pressure would remain 
constant. The flat part of these pressure halts 
lasted from 2 to 5 min, during which time the pressure 
was constant to the accuracy with which the manom- 
eter could be read The first line of table 7 gives 
the observed triple point pressures and temperatures. 


TARLE 7 T ‘ nt temperatures and pressures of the 
j , thanes m the basis cf several analyses ® 
’ , > " y 
4 1 e \/ ! } lig 
minus / a) te i S750 mn 
CHD CE, CHD Cp 
1) 
\/ \ \/ \7 \/ \ \/ 
\! 
t s Hi ‘ ( i Sau ‘ ri 
\ ~ ty van) ‘ a? s SN 
\ j s s ; s si 7 x 3u 
\ j j x 4) 70 r is x 4 8 37 
An } 2. US We 7 70 stil SAS 
*Va of AP mitt ! ( iH 
Cal lated fr i i hl 


2.6. Analysis of the Sample 


Analysis of the deuteromethane samples in the 
mass spectrometer at Columbia University by Pro- 
fessor H.C. Urey and Dr. Marvin Fox revealed that 
the samples were indeed not pure. It was impossible 
to make a reliable analysis at that time because no 
mass spectrometer patterns of pure deuteromethanes 
were available for comparison. However an analysis 
was made based on assumptions regarding the mode 
of formation of ions of various masses. In the inter- 
vening Vears patterns for the deuteromethanes have 
been made in several laboratories and use has been 
made of these studies in determining more accurately 
the compositions of the samples. The various anal- 
yses that have been made are described below. The 
compositions found as the result of each analysis are 
listed in table 8 

Analysis 1. This analvsis, based on patterns made 
at Columbia University by Professor Urey and Dr 
Marvin Fox, was made entirely by a study of the 
abundance ratios of the various masses observed, 
without calibration patterns. The key to. this 
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TABLE 8. Compositions of the deuteromethane samples 


Compouent 


Sample Analysis 
CHy CHD CH,D CHDs CD, 
CHD l 0. 1074 0. 8926 
2 1.00 
3 Loo 
405 1.00 
CH,D | OLS 1917 0. 7968 
2 O76 v24 Trace 
3 12 &Y5 0.004 
4 OW) LOS 785 046 0. 001 
135 S60 OOS 
CHD l OO13 o3oY9 2567 71il 
2 335 652 O13 
; 024 22 283 65S O13 
} O05 1w7 277 AST 024 
5 OS O53 JAS 626 O13 
CP4(2-3 l Ooo! OO44 OSS1 30455 o34y 
2 OSS 023 269 ols 
3 40 OSS 201 ol4 
9.) O37 O48 275 O48 
2-3)8 405 O38 O50 285 26 
(2-4 O3y O53 205 613 
CD, l Os O66 Y284 
3 O20 O68! Vly 
405 OOS O6S w4 
* Values for sample 2-3 a I olate 1 as the mean of samples 2-2 and 2-4 


analysis was the presence in the CD, patterns of mass 
19. Assuming equal sensitivities for the parent 
masses of CHD, and CD,, the CHD,:CD, ratio was 
calculated from the amplitude ratio of mass 19 to 
mass 20. The purity of the deuterium used in the 
synthesis was estimated from this ratio to be 89.26 
percent, and from this purity, the proportions of the 
remaining mixed deuteromethanes were estimated 
in the CD, sample as well as the proportions of all 
the components in the other samples. This method 
was not thought to be very satisfactory because 
slight variations in its application to the various 
sarmples gave different compositions. The compo- 
sitions of the CD, samples thus determined were the 
most reliable, because the ratio of the two principal 
components was obtained with only the single 
assumption of equal ionization sensitivities. 


Analysis 2. The mass patterrs obtained in 1937 
were reimterpreted, using calibration patterps dle- 
termived by Mobler and Dibeler [40] for reference. 
This interpretation wis subject to one error whose 
megritude could rot eccurately be estimated: the 
original pattermrs were obtained on a Nier tvpe mass 
spectrome cer, wl ile the reference pettercs were ob- 
ti ' ed on a Corsolidated Engineering C rporation 
imstrument. Tl e principal vew feature revealed by 
this aralvsis wes 2 much bigher purity for the CH,D 
then had been previously supposed. [> stead of 
conteiming over 10 percent CH, it was found to con- 
taim little or ro CH,. The CH.D, and CHD, com- 
positions were elso rather sh: rply affected; however, 
little corfiderce could he 
the patterns. The amourts of the my jor compo- 
pents of the CD, samples were charged but iittle, 
though more CH,D wes imdiceted to be present in 
semple 2 than at first estimated 


placed in the results for 
these two compourds beceuse of the complexity of 


pe 
m 


of 


po 


pe 


Analysis 3. The study of the deuteromethanes 
made by Schissler, Thompson, and Turkevich [41] 
on a Nier tvpe mass spectrometer became available; 
and this was used also for mterpretation of the origi- 
nal 1937 patterns. This analysis confirmed the high 
purity of the CH,D sample suggested by analysis 2, 
and gave compositions for the CD, samples in reason- 
able agreement with analyses l and 2 However, 
the composition of the CH,D, and CHD, samples 
were not in good agreement with either aralysis 1 
or 2, and this fact suggested the desirability of mak- 
ing a new analysis of the original specimens of meth- 
ane. 

Analysis 4. The original samples were reanalyzed 
using the National Bureau of Standards Consolidated 
EKngineermg Company mass” spectrometer. The 
material used in the previous analysis of the CD, 
sumple 2, fraction from the final purification dis- 
tillation, was longer available. However, the 
neighboring fractions 2 and 4 were analyzed, and 
their compositions are shown in table 8 as samples 
2-2 and 2- The differences are apparently due to 
the slight separation that occurred in the distillation, 
The composition of fraction 3 was taken as the aver- 
age of fractions amd 4. The composition thus 
found does not differ by more than | percent im any 
component from analysis 3. The CH,D, and CHD 
samples showed only a general correpopdence to anny 
of the preceding analyses, differences ranging as high 
as 11 percent. This analysis is subject to two eriti- 
cisms: the analyses were pot xecomparied by cali- 
bration patterrs, and the samples had been standing 
in glass bulbs ¢ losed by stopcocks for approximate sly 
15 vears and might have changed in composition. 
In each of the samples was found more or less nitrogen, 
which presumably entered by leakage of air. Oxygen 


no 


was absent, apparently having been absorbed by 
the stopcock grease The validity of the present 
analyses for the samples as studied in the vapor 


pressure work can be argued from the almost exact 
agreement of analysis 4 of the CD, samples 2-2 and 
-4 with analysis 3, and also from the fact that the 
CD, content of the CHD, did pot show ary increase 


over the vears. The analyses of the CH,D», and 
CHD, samples were still considered to be unsatis- 
factory because the materials used in the previous 


studies of the mass spectra were of rather low purity. 
Analysis 5. Samples of pure CH,D, and CHD, 
prepared by D. H. Rank of Pennsylvania State ( ‘ol- 
lege became available for making lwss spectrometer 
studies in 1953. Beeause these were much better 
samples than had been previously tsed for compari- 
son, apd because it became possible to make a cali- 
bration pattern at the same time an analysis was 
made, the first objection to analysis 4 could be 
removed. <A final analysis of the CH,D, and CHD, 
was made using the Rank samples for comparison. 
The new analyses are closer to analysis 3 than to any 
other, and component differs by more than 4 
percent from that given in analysis 3. Analysis 5 
must be considered to be the most reliable on account 
of the way it was made, though still subject to some 
possibility that the samples have altered in the long 
period since the vapor pressure work was done. 
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2.7. Treatment of the Vapor Pressure Data 


With the information given by the analyses con- 
cerping the compositions of the various samples whose 
properties were measured, it is possible to calculate 
the ee riies of the pure mate ‘rials. In order to do 


this is necessary, in the absence of contrary in- 
te oa to assume the solutions to be ideal In 
this circumstance an observed property, linear in 


the mole fractions of the components, will be given 
in terms of the corresponding properties of the pure 
components and their mole fractions by equations of 
the following type 

AP(1) 


n(l cHpAP n(Densp,AP 


CH3gD 7 CHeDe 


0 
nil cup, Al n(Wen AP CD, 


CHD, 


The number in parenthesis refers to the sample 
under consideration and the subseript compound 
refers to the component, ” is the mole fraction and 
in this instance AY? and A?’® are respectively the 
difference of the vapor pressure of the sample and of 
the pure component from the vapor pressure of CH, 
A total of six simultaneous equations of this type 
were found for the six samples on which analyses 
and vapor pressure studies were made. Using the 
mole fractions given by analysis 5 these equations 


were solved four at a time for the values of A?’° for 
the various components. The solutions are given 
below 
(2) AP CH.D AP(CH,D). 
(b) AP° cup, O.157TAP(CH,D) + 1.16647 
(CH.D,) —0.009AP?(CHD,) 
+ O.001LAP(CD,) 
(c) AP°cup, 0.012AP(CH,D) + 0.5388AaP 
(CH.D.) + 1.602 AP?(CHD,) 
0.023AP(CD,) 
(d) AP°ep, (1) 0.008AP?(CH,D) + 0.040AP 
(CH.D,) —O0.118AP?(CHD,) 
|. OS4AP(CD,). 
(e) AP°cpn,(2--3) 0.043 AP?(CH,D) +O0153AP 
(CH,D,) —0.735AP(CHD,) 
1L.613AP(CD,). 
(f) AP°en,(2—4) 0.045AP?(CH,D) + 0.160AP 


(CH.D,) —0.776AP(CHD,) 
1. 648AP(CD,) 


In applying eq (b) and (c) the CD, sample data 
required for the evaluation of A/’° are those obtained 


from sample 1 which, on account of its greater 
purity, was used in deriving the equations. In eq 
(d, e, f) the CD, sample data to be used is indicated 


in parenthesis. 

It necessary » know the values of the vapor 
pressures of all the pner w at the same temperatures 
in order to use eq (a-f). Advantage was taken of the 
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had been calculated on the basis of the other analyses. 
The short vertical lines in figure 6 indicate the 
spread resulting from the various analyses. a perceptibly better fit, and as it involves an addi- 
: tional constant was abandoned in favor of the two 


Three different functional relations were tried in an 
constant eq (4). 


attempt to fit the data to equations. <A _ plot of 
log!’ pn/ Pu against | T gives a definitely curved line T logywPp/Pa=A— BIT. (4) 
In each case in the liquid range. A plot. of | 

logiol’p Py against 1/7 gave a very n arly straight In the solid range the scatter of the points did not 
line for ach liquid. A calculation of constants for permit a satisfactory choice between the logo?’ p/P 
the equation logyl’p/? y= A— B/T+ CT did not give | plot or the 7 logwlp/Pu against 1/7’; but the latter 
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Fiat RE 6. Ty ) ple points of the deuteromethanes 


was chosen in harmony with the choice for the 
liquid range. 

One fixed point was used in establishing the equa- 
tion for the solid state vapor pressure data. This 
point was determined in such a way as to make the 
discontinuity between the liquid and solid state 
equations for T logy/’?p/P» consistent with the vapor 
pressure equations for solid and liquid CH,. This 
was done by extrapolating eq (2) for liquid CH, to 
the triple point of the deuteromethane and calcu- 
lating at this point log,/’y(liquid) from eq (2) and 
logi/’y(solid) from equation (lb). The term 
Tllogyo?’n (liquid) —logy/’_(solid)] was used to repre- 
sent the discontinuity at the triple point of the 
deuteromethane. It was added to the value of 
T logy’ p/P u(liquid) at that temperature to give the 
value of 7 logy/’p/Px(solid) at the same tempera- 
ture, and this value was used to define the limit of 
the solid state equation. It may be noted at this 
point that the revision of the denominator constant 
of the Antoine equation for solid CH, has a significant 
effect in this calculation, as shown in the following 


tabulation 


Tllogio Px (liquid) — logy) Py (solid) 


Source of function 


CH.D CH.D, CHD CD, 

Equation (la) 0.157 0.291 0.407 0.500 

Equation (lb) (re- 145 . 280 . 396 . 489 
vised). 

Estimated from data 147 . 276 . 392 . 487 


In every case the revised equation gives a value 
for the discontinuity that is more nearly consistent 
with the experimental points than does the unrevised 
equation. 


| 


| 


Table 10 gives a complete list of the constants for 
eq (4) fitting the data calculated on the basis of the 
various analyses. In the last column e represents 
the mean deviation in millimeters of mercury of the 
data from the equations. The deviations were 
calculated for equations resulting from analysis 1. 
With the exception of the CD,, they should be 
approximately the same for the other analyses. 
The mean deviation, +0.027 mm, of the CD, data 
shown for analysis 1 is the deviation of the data of all 
the CD, samples from the equation fitted to the 
average of all. However, much of this represents 
systematic differences between the samples, as 
sample 1 data lie consistently below the sample 2 
data. In the final treatment of the data based on 
analysis 5 the equation was derived to fit only the 
data from sample 1. For this sample the mean 
deviation was 0.009 mm. The constants listed in 
table 10 show clearly the differences resulting from 
the different compositions attributed to the samples. 


TABLE 10. Constants for the equation T log Pp/ Pu A B/T 
Solid range Liquid range 
Analvsis Methane e (mm 
is Bs \1 Br 
CHD 1.424 124.7 1481 1445 40.015 
, CH,! 2.842 234.2 2781 2546! +.017 
CHI $152 320.3 3.773 328.5; +.018 
CD, (all 5.13 74.8 1360 419.1 +. 027 
CH,D 
; CH)D 
= } CHI 
CD, 
CHD 1.324 129.2 
3 CH,! 2. 632 242. 6 
; | CHD 3. 942 343.4 
| CD, (all 5. 206 425.1 
| 
CHAD 1.259 110.2) 1.3275 | 120.5 
, | CHI 239) 195.7 2.742 | 253.5 
| CHI $576 361.7 4.101 355 
CD, (all $915 355.3 5.138 | 419 
ch. 5.357 | 396 
CH, 1. 260 110.2 1. 328 129. 5 
CHD 604 222.2 2.671 | 245.4 
| CHI 41.452 351.7 3.969 343.8 
CD, 5.529 410.5 5.159 | 421.1 +. 009 
i 
| A further check on the internal consistency of the 
| data was made in the following way. Using the triple 


point pressures and temperatures of the four deuter- 
omethanes shown in table 7, the function 7 logioPp/ 
P,, for the solid at this pressure was calculated and 
| compared with the value taken as a fixed point in the 
| derivation of the solid deuteromethane equations. 
The values found are compared below. 


Sample 


i logioPp Pu 


From equation 
for solid 
deuteromethane 


From triple 
point data 


CD, 0. 95 0. 958 
CHD 54 543 
CH,D, 24 . 230 
CH,;D . 07 . 041 


50 


wel 
anc 


e 
le 


r- 


D 
id 


1e 


The agreement is good. Only in the case of CH;D 
is an appreciable error observable. If the triple 
point of CH;D were 84.45 mm instead of 84.52 mm 
the values for T logyPp/P H would be in agreement. 
This difference in triple point pressure could be ac- 
counted for in either of two ways, (1) there may be a 
2 percent impurity of CH, in the CH;D sample, 
which is not an unlikely possibility, or (2) an error 
of 0.07 mm was made in reading the triple point 
pressure 


2.8. Discussion of the Vapor Pressure 


The vapor pressure ratios Pp/P?y show certain 
regularities of behavior which are worthy of mention. 
The substitution of each successive protium atom by 
a deuterium atom causes a nearly constant change in 
the constants A and B. A slight but consistent 
deviation from a linear increase of A and B occurs 
as may be seen by reference to figure 7. It has been 
postulated as a working approximation in estimating 
the vapor pressure of the middle member of a series 
of increasingly deuterium substituted compounds, 
such as the waters, that the vapor pressure of the 
middle member (HDQO) is a geometric mean of the 
vapor pressures of the extreme members (H,O and 
D.O). Within the limits of experimental error the 
few determinations on HDO support this hypothesis, 
and a theoretical justification has been presented 
[42]. In the case of the methanes a similar hypothesis 
would require that the ratios Pp/P?y should form a 
geometric progression. This would require A and B 
to form a linear series. The nonlinearity of the 
constant A indicates that deviations from a geo- 
metric progression do exist though they are small. 
The deviations from a geometric progression are also 
evident in the fact that the ratios Pp/Py are not 
unity at exactly the same temperatures; but the 
temperatures at which this ratio is unity are suc- 
cessively lower the more highly deuterium substituted 
is the methane. The melting points also form a 
nearly linear series as may be seen in figure 6. The 
deviation from linearity appears to be real, though 
slight. 

Because of the form of the equations representing 
the vapor pressure ratios, it appears that from the 
temperatures at which measurements were carried 
out up to the critical point the vapor pressure of the 
deuterium compound will always be greater. If the 
same form of vapor pressure relation is applicable 
throughout the liquid range then it appears that the 
vapor pressure ratios will have maxima at or below 
the critical points. The maximum vapor pressure 
ratios indicated by the equations range from 1.009 
in the case of CH,D—CH, to 1.037 in the case of 
CD,—CH,. 


2.9. Ideality of the Solutions of Isotopic Isomers of 
Methane 


In view of the complexity of the two mixtures which 
were used as CD, samples in the work of Brickwedde 
and Scott, and the difficulty of obtaining an accurate 
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NUMBER OF DEUTERIUM ATOMS PER MOLECULE 
FiGgureE 7. Constants of the equations T logy Pp/Pa= A—B/T. 
analysis, it appears that only the roughest estimates 
can be made of the magnitude of deviations from 
ideality which they exhibit. Two estimates have 
been made. 

One estimate has been made of the deviations 
from ideality to be expected if the deuteromethanes 
form regular solutions, according to Hildebrand’s 
criterion, that is, have ideal entropy of mixing. For 
this type of solution the activity of a component 1 
is given by eq (7). 


RT log, a,=RT log, X,+ Vi¢3[(AM,/V,)'” 
(AH,/V2)'/7] (7) 


is the latent heat of vaporization per 
cubic centimeter, @ is the volume fraction and V 
the molar volume of a pure component. Using the 
following experimental values the activity coefficient 
of CH,D in its 10 percent solution in CH, was calcu- 
lated to be 1.000058. This calculation places the 
possibility of finding deviations from ideality of this 
type beyond experimental precision. 


where AH/V 


Sample AH{21] V [22] 
cal cm 
CH,. . 2, 036 33. 63 
CH;D_- 2, 050 32. 4 


Another extremely rough estimate of the deviations 
from ideality in the solutions rich in CD, has been 
made by comparison of the vapor pressure of CD, 








calculated from samples 1 and 2-4. 
that sample 2-4 gives values always higher than 
sample 1. This suggests that sample 2, which con- 
tained more admixed CHD, had a vapor pressure 
higher than it would have had if the solution were 
ideal. If we consider the vapor pressure of pure 
(‘D, to be that calculated from the purer sample 1, 
then y, the activity coefficient of CD, in sample 2, is 
given by Pep, (sample 2-4 calculation) =yPcp, (sam- 


ple 1 calculation). The value of y thus calculated 


averages about 1.00028 in the liquid range, and does | 


not vary much with temperature. In the solid range 
it increases from 1.0008 near the triple point to 
1.0035 at 77°. ‘These values appear to represent the 


activity coefficient of CD, in a solution of about 70 


percent CD, as compared with its activity coefficient 
in a solution of about 92 percent CD,, the assumption 
being made that the vapor pressure measurements 
are not systematically wrong for either of these 
samples. 

The relative volatilities of the deuteromethanes 
with respect to CH, have been calculated at the nor- 
mal boiling point of CH,; they are: CH,D, 1.0035; 
CH,.D,, 1.010; CHDs,, 1.019; CD,, 1.029. 
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Data on the Atomic Form Factor: Computation and Survey 


Ann T. Nelms and Irwin Oppenheim 


This paper presents the results of calculations of atomic form factors, based on tables of 


electron charge distributions computed from Hartree wave functions, 
Computations of the form factors for five elements 


atomic numbers, 
arsenic, and mercury 
cated, 
is reviewed, 
sparse experimental data are made, 


1. Introduction 


The atomic form factor is of interest in the caleu- 
lation of Rayleigh scattering of radiation and coher- 
ent scattering of charged particles from atoms in the 
regions where relativistic effects can be neglected. 
The coherent scattering of radiation from an atom 
consists of Rayleigh scattering from the electrons, 
resonant electron scattering, nuclear scattering, and 
Delbriick scattering. When the frequency of the 
incident photon approaches a resonant frequency of 
the atom, large regions of anomalous dispersion occur 
in which the form factor calculations are not sufficient 
to describe the coherent scattering. 

The cross section for Rayleigh scattering from the 
bound electrons of an atom with a deflection ¢,' per 


steradian dQ, is 

/ ro 2a) lf \/2 10 

da=~ (1+ cos" @) |f(qg)\° dQ, (1) 
where r=—@/me?—2.82“10-" em is the classical 


electron radius, f(g) is the atomic form factor, q is 


the change of momentum of the photon, and 
lq (hv/e)2 sin(@/2), where y is the frequency of the 


incident photon. Terms of order v/c, where v is the 
initial velocity of the bound electron, are neglected 
in the derivation of (1 

The cross section for the coherent scattering of 
charged particles with a deflection ¢, per steradian 
dQ, is, in Born approximation: 


do [Z f(q)] cosec” ¢/2 dQ, 


ef )) 
(2 
2 Mi" ; 
where ¢ is the charge of the scattered particle, M is 
its mass, and ¢ its velocity. Z is the atomic number 
of the atom from which the particle is scattered. 
Equation (2) is a good approximation for all angles 
if the kinetic energy of the incident particle is large 
and Z is not very large. For lower kinetic energies, 
eq (2) is valid for small-angle scattering. 
more usual notation ¢=26, where @ is the Bragg angle. 


| 
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for a wide range of 
carbon, oxygen, iron, 


are presented and a method of interpolation for other atoms is indi- 
\ survey of previous results is given and the relativistic theory of Rayleigh scattering 
Comparisons of the present results with previous computations and with some 


The form factor for an atom of atomic number Z 


is defined as the matrix element: 


¢ ig: r/h 0 ) 


where 0 denotes the ground state of the 


> Z 
fg=(0 >> (3) 


l 


| 


atom, and 


r, is the vector distance of the jth electron from the 
nucleus. For a spherically symmetric atom 
: oF sin kr 
f(q) p(r) dr, (4) 
‘ 0 kr 
where k=47(v/c) sin (¢/2)=q/h 


and p(r) is the charge distribution of electrons. 

The basic data for the form-factor calculations of 
this paper were the effective charge distributions in 
atoms recently tabulated by Freeman and others. 
The effective charge, Z,(r), may be defined in terms 
of the electron charge distribution 


PZ,(r) 


(1/477) dy? 


p(r)- (5) 


Because of the oscillatory nature of the integrand in 
eq (4), the calculation of the form factor must be 
done in a careful manner. 


2. Historical Survey 


The initial attempt to evaluate p(r) numerically 
was made by Hartree [1],’ using the Bohr model of 
the atom. Atomic fields for a few elements were 
obtained from a numerical analysis of optical and 
X-ray spectra from which core orbits were calcu- 
lated. The results were generalized so that the form 
factor for an arbitrary ion could be calculated. 
Important data for the analysis of crystal structures 
were obtained in this manner, but the calculated 
form factors were, in general, larger than the ob- 
served values. 


2 Figures in brackets indicate the literature references at the end of this paper, 








The advent of the wave mechanical theory of 
scattering by Wentzel and Waller [2] confirmed the 
general formulations (1) and (2) of the scattering 
process. Charge distributions were now to be de- 
termined from solutions of the appropriate Schréd- 
Inger equations An exact solution is obtainable for 
the hydrogen atom but not for many-electron atoms. 
For a one-electron system 


* l ; 6 


pl) far-y* (rw 
where ¥(7) is the Schrédinger wave function. Some 
of the approximations used for many-electron atoms 
W ill be discussed below. 


2.1. Pauling and Sherman Model 


Pauling and Sherman [3] obtained form factors 
for all atoms and some ions from hydrogen-like 
eigenfunctions with sereening constants character- 
istic of each subshell. Very good results for small- 
angle scattering in all materials and for inter- 
obtained. 


mediate angles in low-Z materials are 
The method tends to overemphasize the shell 
structure of atoms, and one would expect: more 


accurate results to lie between the Pauling-Sherman 
curves and the Fermi-Thomas curves 


2.2. Hartree Self-Consistent Field Model 


The most accurate extensive computations of 
wave functions of many-electron atoms are based 
on Hartree’s [4] self-consistent field method. This 
is an independent particle model in which each 
electron is assumed to be in the field of the nucleus 
and an average field due to the other electrons. 
Thus the charge distribution can be considered as 
being a sum of one-electron charge distributions 


Z Z 


pl >> oj) =4ar? SO VW 0 7 
where vy; r) is the wave function for the jth electron. 


This treatment is applicable to ions, as well as 
atoms. Fock [5] has generalized the Hartree scheme 
by including the effects of exchange, but few com- 
putations of Hartree-Fock wave functions are 
available. Slater [6] has shown that the Hartree 
and Hartree-Fock wave functions arise out of 
variational treatments in which the wave function of 
the many-electron atom is assumed to be a product of 
individual electron wave functions. 


James and Brindley [7] have utilized Hartree 
wave functions in the calculation of atomic form 
factors. They tabulate results for 1<Z<37 and 


O<(sin @/2)/A<1.1, using an interpolation method 
for elements not given by Hartree. Their results 
are presented in tables of form factors for subshells 
of electrons versus [1/(Z—s)]|(sin $/2)/A], where 
s is a screening parameter determined empirically 
as a function of Z and the subshell of electrons 
considered. The atomic form factor is obtained by 
summing the form factors for the various groups of 


electrons. The method provides a suitable means of 
interpolation. The calculation of the form factor 
involves the numerical integration of a product of an 
oscillating function and the charge distribution for 
each subshell. Computational errors are inherent 
in this procedure and when the results are summed 
to give atomic form factors appreciable errors may 


occur. 


2.3. Fermi-Thomas Model 


Fermi [8] and Thomas [9] treated the electronic 
cloud of an atom by Fermi-Dirac statistics and ob- 
tained a spherically symmetric electron atmosphere 
having a density 


Z o(.\) 7 
P 7a | P 9 ) ' ») 
where uw=0.47X107°°*& Z em, V=r/p, and 
o(.V)/N=(u/Ze) V(r), where 
Vn=“ ef alr) py 
/ , 


is the potential at 7 due to the nuclear charge Z and 
the electron charge distribution p(r) 


Form-factor calculations have been made as a 
function of the universal variable 
u—4aru sin(d 2)/Xr, 


where \ is the wavelength of the incident radiation. 
The calculations have been applied to neutral atoms 
and positive ions and give smooth electron distribu- 
tions that average the effects from each electron. 
This approximation is fairly good for high-Z materials 
over most of the electron distribution but leads to 
large discrepancies with experiment at small and 
large angle scattering. The inherent limitations of 
the statistical method are the following: 


(a) The effects of atomic shell structure are 
smoothed out. 

(b) A poor distribution is obtained where the 
potential is small (at large distances from the 


nucleus) and where it changes rapidly (close to the 
nucleus). 

A discussion of the inadequacies of this model for 
incoherent scattering is given by White [10]. 

Gaspar [11] has noted that reduced atomic form 
factors determined by the self-consistent field method, 
can be regarded as almost universal functions of the 
Fermi-Thomas variable 


0.8853 


Z l 


u ak=0.8853 L', (9) 


where k=4x(sin $/2)/A=q/h, and ap is the Bohr 
radius. Asis evident from figures 1 and 2, this is true 
only for values of (’<3. There are appreciable 
differences at higher values of (° due to the shell 
structure of the atoms. 


*/Z 


= t+ -- 
see figure 
cons Suma sad 
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3. Relativistic Effects 


Relativistic effects are of importance for heavy 
atoms where the inner electrons have relativistic 
velocities and for scatterings in which the momen- 
tum change of the photon is of the order of or greater 
than Calculations in this report have been 


carried up to momentum changes less than or equal 
) 


mic. 


to ite 

Franz [12] has shown that the relativistic expres- 
sion for scattering of a high-energy photon by a bound 
electron reduces to the form-factor calculation when 
q<me and the velocity of the electron is non- 
relativistic, i. aZ/n 1, where @ is the 
fine-structure constant, and 7» is the principal quan- 
tum number of the electron. His method of cal- 
culation utilizes plane waves for the electron wave 
functions in the intermediate (The Dirac 
perturbation theory describes scattering as a two- 
step process, absorption followed by reemission or 
emission followed by absorption; thus one must 
consider the state of the electron between the two 


LS 


e., (v/c)*~ 


states. 


steps.) The effect of electron binding in the inter- 
mediate state is thus neglected. For his computation 
of the form factor, Franz uses a Fermi-Thomas 
distribution. 

Note that F iat 2). whiel valid for seat ng of high ener 
phot nly, v a factor of 2 due to a mistake in transforming variabl 
of int atic ver ingl f i ! Thus we have for the total « ectio 
eot “Fes . 
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Levinger [13] has considered corrections to the 
form factor calculation in the near relativistic region 
but has neglected the effect of binding in the inter- 
mediate states. Bethe [14] has confirmed Franz’s 
result and has shown that for values of q larger than 
the characteristic momentum of a A’ electron 
meZ/137, Dirae wave functions should be used for the 
computation for A electrons. Rohrlich and Rosen- 
zweig [15] have extended Bethe’s results to a treat- 
ment of L electrons. 

Brown and Woodward [16] have investigated the 
effect of binding in the intermediate states and find 
further corrections to the form-factor calculations 
of the same order of magnitude as Levinger’s. An 
extensive project by Brown et al at Birmingham in 
the exact calculation of Rayleigh cross sections for 
selected atoms and energies is now under way. 
Results for the scattering of 0.32 me? gamma rays 
from the A electrons of mercury [17] indicate that 
the form factor results are a reasonably valid descrip- 
tion of the dispersive scattering but not of the total 
scattering which includes contributions from absorp- 
tive scattering. Under the conditions of the calcula- 
tion, the absorptive scattering is approximately \% 
of the dispersive contribution. 


4. Form-Factor Calculations 


The form-factor calculations of this report are 
based on the published numerical data of Freeman 
[18] and others [19], who calculated the effective 
charge distribution in atoms, Z,(r), from the Hartree 
self-consistent field model.4 Hartree-Fock data were 
not used because they were unavailable in quantity. 
It is well to note that, although the radial charge 
distributions of the Hartree and Hartree-Fock 
models differ considerably, the corresponding form 
factors may not. Brindley and Ridley [20] have 
computed form factors for K*, Cl-, and Cu*, using 
both Hartree and Hartree-Fock wave functions. 
Their computations, which cover the range from 
k=0 to k=2rn, indicate that the form factors with 
exchange are higher than those without exchange. 
The maximum difference in the form factors is of 
the order of 3 percent in this small-angle region. 
A trial calculation in the present work for Cu* at 
k=10 indicates a difference in the form factors of 
0.6 percent between the Hartree and Hartree-Fock 
results. 

The effective charge distribution Z,(r 
by the relation (5 


is defined 


PZ, \T) 
dr? 


(ser) 
rr 

, tar 
where p(r) is the charge density of electrons in the 
atom. Equation (4) may be rewritten, by means of 
partial integration to give 


a) 
’ P . , 
f (I) Z— | Z,(r) sin (kr) d(kr), (10) 
0 
An upper limit of 70 (q=me/2) was established for t calculation that 
relativistic effects could be neglected wit! t serious er 


where kK=q/h. In the tables of [18], r is given in | where the sum over m denotes a sum over the various 
units of the Bohr radius, and & must be expressed in | intervals for which Z,(r) was fitted, ro=0 and 


the reciprocal units Tinmay = ©, and the superscript (m) on the coefficients 
a; and a@ denote the values for the appropriate 
kag= (42/X) sin (¢/2)ao=6.635 sin (¢/2)/X. intervals. An analytic integration was obtained 


in this way for 0.5< k < 20 (fig. 3 and table 2). 

The numerical data. Z.(r). for each atom were | The maximum error introduced by the fitting 

fitted empirically by an expansion of the form | procedure in the calculation can be easily estimated 

by assuming that the deviation of the fit from the 

ne tabulated value is in phase with the oscillating 

e~or Doar’. (11) | sin(Ar). An error of 3 percent is estimated in this 

way. However, it is improbable that the actual 

error is of this magnitude. A more realistic value 

would be of the order of 1 percent. Although this 

calculation procedure is straightforward, the length 

of the calculation and inability to make general 
checks make it a tedious operation.° 


Functions of this form were chosen so that direct 
comparison could be made with the effective charge 
of hydrogen type atoms 


Z,(r)=e-"(Z+rZ). (12) 
rh —e —e : TABLE 2. Form factor, f Z, tabulated from eq | 13) 
The fitting procedure for nonequidistant values of 7 
is described in detail in [21]. The Z,(r) data were 


fitted in several intervals separately. This was done “ 
to obtain a higher degree of accuracy without using ‘ - - o- 
expansions involving large numbers of terms. In | 
this way, six-term expansions (n=5), at most, were | 0.0 1.0 1.0 1.0 1.0 1.0 
necessary to fit the data. Values of the coefficients | e “— = al “a. = 
for each interval for the five atoms considered—C, | 2.0 415 . 508 . 5R2 675 700 
O, Fe, As, Hg—are listed in table 1. It is, perhaps, | 4.0 ‘Ont “240 * 30 “436 “530 
only in the first interval that one can make a signifi- | 2 = ‘oo. - _ 
cant comparison with the effective charge of hydro- | 8.0 . 160 . 150 217 21 — 
ah . " 10.0 104 . 118 . 182 . 178 250) 
gen-type atoms. The aecuracy of the tabulated | 20.0 . 0217 0386 . O785 . 0890 .109 


Z,(r) and of the fitting is noted also in table 1. In 
all cases effort was expended to make the fit as accu- 
rate as the data. The behavior of the form factor for small k, k- 

0.5, was determined by expanding sin sr in powers 


of kr: 


Substitution of (11) into (10) results in 


f (I) DOS, cd 
y =1—(k Z) 22 1) | fh 242 = (2 a Z,(r)r°"~ "dr. (14) 


oO! (''m ” | 
- ae ° » 
"iis sin (kr)dr, (1: , 
i ) Note that a graphical procedure would not have been suitable in this range 
. 1 because of the oscillating nature of the integrand 


raBLe 1. Coefficients in the expansion Z"(r)=e7-@"™ * Yar 
i 


Maxi Maxi 
mum mum 
7 Interval, a t¢ ay a? 13 a4 ag error of error in 
calcula reference 
tions [18] 
{0 to 0.1 2 47 6 0.01 | 
Cc 1 to .72 2 6.04511 + US448 18. 60504 11. 06320 Ol 0. 01 
| .725 to 1.§ 4. 84325 1. W406 1. 23017 0. 158 o | 
{0 to 0.1 2.77 s 2 | 
oO = ] to .8 2 7. 673 $. H22 13. 42 10. 34 Ol » 0.02 
is ¢t 1.9] 7. 764 0.917 1. 884 0. 395 OOS | 
} 
(0 to 0. OR1 § 7283 Pai — 20). 3322 256). O535 61. 5426 02 
} .OSl to .73 2. 35806 24. 68559 —15. 43486 62. 94042 548. 52436 1032. 86606 610. 402954) ll 
Pe ) 735 to 2.0 1. 3966 24. 11512 —24. 18032 14. 23867 2. 27092 0001 9.99 
120 to 1. 07624 5 1507 6 1, 28879 —(), 07209 0. 00950 003 | 
10 to 0.1 { 33 20. 55955 671. 12959 03. H4918 Ol 
ae ian to 7 2. ¢ 34. 28291 SY. 423502 464. 62137 1124. 45247 1202. 72800) 551. 52777 Os 0.03 | 
7 01.8 1.27713 39. 90909 —94. 09721 122. 21415 79. 41765 26. 37782 3. O417 O02 . | 
L8 t 1. 14162 9, 79098 0. 14061 0. 24590 (0. 038708 003 
0 to 0. 106 5. 14277 SO) — 59. 11228 2524. 22231 124272. 83758 146263. 92891 6462387. 83717 03 | 
O6 to .478 $ 91425 45. 99829 637. 7033 6844. 959047 31986. 00973 — 67636. 26526 52050. 76587 Ol | | 
He 178 to 1.96 1. 68832 14. 19283 24209358 | —469. 93545 420. 91024 178. 44611 29. 20347 04 p Gee | 
li 0 te 1. 09654 31. 50632 ~—23. 61074 ¥. 58349 1. 72085 0. 13744 — 0. OO386 Ol 
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FiGuRE 3 Form facto n the ange 0 I< k« “ 
At the origin, one obtains: 1.0 
~~ 
} ‘ ——— — — 
; dfik) tS 
(0 Z and ( ) 0. Y 
dk SN 
Near the origin: | = SS ie 
. )) j 
f(k Z+k | Z,(ryrdr. D | | | ; = ee | Hg 
. i oll 
Nosh a a 
The integral in eq (15) was evaluated graphically 0 
with the tabulated Z,(r) (see fig. 4). KN Yas 
It is of interest to note the connection of the form Ss = i\ : 
factor for small-angle scattering with the expression } Fe 
for the gram-atomic diamagnetic susceptibilities. — - ~ 
The expression for the susceptibility, Xdia, 1S Yc 
| — 
“= 8X 107*)ao72(0! S *\r,!7/0). 16 90L_. ie 
; a ast 0 , 2 3 4 5 
k=47(0,/r)SiN(9/2) 
where ad) is the Bohr radius. For spherically 
. Figure 4 Small-angle approrimation of the form factor per 
symmetric atoms, we have 
‘ efectron, 
0) S3ir,;!7'0) | Z(ryrdi \7 ’ 
‘ Z,(r) are for free atoms, and since the experimental 
and from (15 Xa... as well as most data on the form factor from 
je? X-ray and electron small-angle diffraction, is meas- 
f(k Z 2 6 (1.25 10 Xia (1p. lS ured for atoms either bound in molecules or crystals, 
) . ° 7 . 
a comparison between the computation and experi- 


ment is not quite direct. Brindley and Hoare [22 
have discussed the difficulty in evaluating sus- 
ceptibilities from experimental data. Nevertheless, 
the diamagnetic susceptibilities computed from the 


The integral in (17) involves a precise knowledge of 
the effective charge, or the wave function at large 
radii and is therefore strongly dependent on the 
chemical binding of the atom. Since the tabulated 
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values of Z, used in this paper agree with the experi- 
mental values obtained for the atom in chemical 


combination to within 3 to 10 percent (see table 3). 


TABLI 3. Gram-alt ride susce plibilit /, X4 10' 
7 ( ul I erimental 
Cc ‘ 
2) OS j 
Fe 
As 2») 2% 
Hy } is. 1 
E. C. Stoner. M I ( Ltd Essex 8 t, W. ¢ 
Handbook of ¢ ( Rubber Put 
( Cle land, 
For > 20, rough calculations were made because 


of the uncertainty of errors due to relativistic effects. 
For large & a useful expansion of (10) is 


f(h: | G(r) cos (kr) di 


t\"G™() 
! > & ‘OS [A ( 
: (7). eos (kr) dh 
1 G(r) cos (kr)dr- - | G(r) eos (kr)dr 


, = } G(r) cos (kr)dr 


| : t\*G'™ (r) 
r% >> ‘os (kr)dr (19) 
SJ. CE) am cos 


TN\k 


100 





where 
] dZ,(r) 
Z dr 


dl" 


G(r) = and G\” (r)=.~ Gir). 
dr" 


The coefficient of the G° 
approximation used was 


rk) a | rik 


e 


integral is zero. The 


Gir) cos (kr) dy T " G(r) cos hy dr 


: (7) l" | G r) cos (hy dr. 20) 


~ 4ylhE 


As is to be expected, only regions close to the nucleus 
contribute to the form factor for large momentum 
changes. These computations were compared with 
a hydrogen-like approximation, using Slater sereen- 
ing constants for low-Z materials. For high Z a 
comparison with a Thomas-Fermi calculation was 
made. The comparisons were informative for gen- 
eral trends only and could not be used as a test of 
accuracy. 

An interpolation scheme is given in figure 5 over 
the range for which calculations were made, 
0.5<k<20. Tabulations of the form factor by 
James and Brindley, which are based on Hartree 
wave functions, were used as a guide in regions where 
we had little calculated data and as a general test 
of interpolation in other areas. The curves are 
plotted with 1l-percent accuracy, and interpolation 
can vield results within 2 percent. A more direct 
test of the interpolation was made by comparing 
the form factor of Cu*(4=10), obtained graphically 
from published data [23] of Z,(r), (f=5.26), with 


the interpolated value for Cu (f= 5.23). 


10 100 


f 


Figure 5. 


Form factor as a function of Z in the range O<k < 20, k 


tm(do/dr) sin (b/2). 
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5. Comparisons With Other Calculations 


Com parisons of the present results with the numer- 
ical results of the approximations discussed in section 
3 are given in figure 6 for C, O, Fe, As, and He. 
Figure 7 shows a comparison of the results for Sn and 
Al obtained from the interpolation curves, figure 5, 
with the calculations of section 3. 

As is to be expected, the data of Pauling and 
Shern an agree best for low Z with the present more 
exact calculations. In general, the oscillations in 
the Pauling-Sherm an form-factor curves due to shell 
structure are larger than in the more exact 
calculations. 

The caleulated curves oscillate about the Fermi- 
Thomas curves, which average out the shell structure. 

Comparison with the results of James and Brindley 
shows agreement over most of the range. The 
present form-factor computations are higher than the 
James and Brindley calculations for small and 
moderately large values of /, especially for C in the 
' range 0.655<k<2. The results reported here show 
Interpolated close agreement with the extensive computations by 
MeWeeny [24] and others on the form factor of C. 

The present authors judge the results of this paper 
to be more accurate than those of James and Brindley 


Pauling and Sherman 


James and Brindley 





S 1 10 
k= 41(a,/a)sin (6/2) results jrom other calculations. 


FIGURE 7. Comparison of interpolated form factors with 
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because the present computations utilize more recent 
charge-distribution data [18] and also avoid the 
errors inherent in the sum over subshell contributions. 

The ratio of the cross sections fer coherent and 
incoherent scattering should be a smooth function of 
Z. The value of this ratio for C obtained from 
MeWeeny’s calculations and from our own fits in 
well with the values for other light elements whereas 
the values obtained from James and Brindley’s 
calculations do not. This anomaly of Janes and 
Brindley’s data for C was apparent even in an anal- 
ysis of total absorption coefficients (scattering and 
photoelectric) and contributed to a stimulation of the 
present study. 

For high Z materials, Z>26, comparisons are 
made with the Fermi-Thomas and Pauling and 
Sherman calculations. As was expected, the more 
exact calculations approach the Thomas-Fermi 
results as Z increases. Even for Hg, however, some 
effect of shell structure is apparent. A recent com- 
putation of the form factors for Hg** by Henry [25] 
involving a numerical integration of Hartree wave 
functions is plotted for comparison in figure 6. 

In figure 7, interpolated form factors are compare: 
with those of Pauling and Sherman for Sn and with 
those of Pauling and Sherman and of James and 
Brindley for Al. There are deviations of approxi- 
mately 17 percent from the data of Pauling and 
Sherman for k= 2.64. 


6. Comparison With Experiment 


In the range for which our results are valid, i. e., 
nonrelativistic momentum changes, and in’ which 
there are significant differences in the results of the 
various form factor calculation, (> 3, there are very 
few experimental data with which the calculations 
may be compared. The early X-ray data give 
information only for small momentum changes, 
(<3 in figures 1 and 2, where the Hartree and 
Pauling-Sherman approximations approach the “uni- 
versal”’ Fermi-Thomas form factor calculation very 
closely. Recent experiments by Wilson [26] and 
Mann [27] were not used because the momentum 
changes were larger than me/2, the limit established 
for the present calculation. 

Comparisons are made with the experiments of 
Storruste [28] with 0.411-Mev incident photons on 
Pb and Cu (fig. 8). The experimental error is 
indicated by the size of the experimental points. 
Satisfactory agreement is obtained throughout the 
range though the experimental points for momentum 
changes greater than me 6 lie above the theoretical 
curve for Pb. 

The form factors that have been computed here 
are for Rayleigh scattering from the electrons of free, 
spherically svmmetric atoms. MeWeeny [24] has 
made extensive studies of the effects of asymmetry 
on the form factors, particularly for C. Using 
analytic wave functions of the Duneanson and 
Coulson type [29], he finds effects of the order of 3 
percent on low- and medium-angle-scattering factors. 
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form factors with experimental data, 


Experimental determinations of coherent scatter- 
ing are almost always made with atoms bound either 
in molecules or in solids. Interatomic binding will 
distort the electron cloud producing a modification 
in electron charge distribution and therefore in the 
form factor. This binding will affect the small 
angle scattering of photons whose wavelengths are 
of the order of atomic dimensions, (1 to 10 A). For 
higher energy photon scattering only the electron 
distribution near the nucleus, which is hardly in- 
fluenced by interatomic binding, is of importance. 
Brill [80] has discussed the effect of interatomic 
binding, particularly in crystals, and finds effects of 
the order of 4 to 5 percent for [sin (@/2)]/A~0.5. 
MeWeeny has found no significant difference between 
form factors computed for the valence states of free 
atoms and the form factors for these atoms bound in 
crystals. 


Furthermore, as mentioned above, coherent scat- 
tering from atoms consists of Rayleigh scattering, 
resonant electron scattering, Delbriick scattering, 
and nuclear scattering. Resonant electron scatter- 
ing produces large regions of anomalous dispersion 
when the frequency of the incident photons ap- 
proaches a resonant frequency of the atom. Parratt 
[31] has shown that within a large frequency range, 
near and below the K edge, e. g., 0.1 AC A< 100 A for 





Cu, there is practically no region of normal disper- 
The form-factor calculation is not sufficient to 
obtain the coherent scattering in this region. 

For a comparison with experimental results at 
higher frequencies, the cross sections for the Ray- 
leigh, Delbriick, and nuclear must be 
known, as well as the phase relationships between 
the scatterings by each Rayleigh and 
nuclear scattering will constructively interfere up to 
several million electron volts while Delbriick scatter- 
ing is out of phase with them in this range. 

The inapplicability of the Born approximation (2) 
to the seattering of electrons by heavy atoms has 
been considered by Glauber et al. [32]. Tabulations 
of form factors corrected for deviations from the 
Born approximation have been given by Ibers and 
Hoerni [33] 


sion 


processes 


process. 
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